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PREFACE 


This  final  report  covers  the  work  accomplished  during  the  period  from 
March  1,  1994  to  February  28,  1997  under  the  Air  Force  Office  of  Scientific 
Research  (AFOSR)  contract  F49620-94-0137.  The  AFOSR  technical  monitor  for 
this  contract  was  Dr.  Charles  H.  Ward. 

Gamma  titanium  aluminide  alloys  are  potential  structural  materials  in  the 
aerospace,  power  generation  and  automotive  industry  due  to  its  superior 
mechanical  properties  over  the  conventional  superalloys.  However,  wide 
spread  application  of  these  materials  necessitate  improvement  of  their 
overall  mechanical  properties  and  reduction  of  their  cost  by 
optimizing/exploring  processing  methods.  To  that  effect,  spray  atomization 
and  deposition  of  gamma  titanium  aluminide  alloys  were  investigated  in  the 
present  study,  with  the  following  specific  objectives.  First,  to  determine  the 
feasibility  of  using  a  spray  atomization  and  deposition  approach  to  S5mthesize 
y-TiAl  with  improved  hot  workability.  Second,  to  develop  a  fundamental 
understanding  of  the  influence  of  spray  atomization  and  deposition 
processing  on  the  microstructural  characteristics  (e.g.,  grain  morphology, 
segregation,  porosity)  of  y-TiAl.  Third,  to  elucidate  the  influence  of  these 
microstructural  characteristics  on  the  relevant  deformation  and  fracture 
mechanisms  that  govern  the  behavior  of  y-TiAl.  In  order  to  achieve  the 
above  mentioned  objectives,  both  experimental  and  numerical  methods  were 
employed.  Experimental  approaches  included: 

1) .  spray  atomization  and  deposition  of  y-TiAl  alloys,  the  prerequisite 

to  the  success  of  this  program; 

2) .  mechanical  testing,  such  as  tensile  testing,  compressive  creep 

testing,  and  tensile  creep  testing  on  spray  formed  y-TiAl  alloys;  and 

3) .  characterization  of  microstructure  and  phase  compositions  of  as- 

received,  spray  formed,  and  crept  y-TiAl  alloys,  using  optical 
microscope,  image  analysis,  density  measurement,  SEM,  TEM,  and 
X-ray  diffraction  analysis. 

These  three  closely  related  aspects  were  aimed  to  experimentally  establish  the 
relationship  between  processing,  microstructure  and  mechanical  properties. 
Numerical  analyses  include: 


1) .  Modeling  of  momentum  and  thermal  behavior  during  spray 

atomization  of  y-TiAl; 

2) .  Investigation  of  thermal  residual  stresses  in  spray  formed  y-TiAl 

using  finite  element  method;  and 

3) .  Study  on  creep  and  boundary  sliding  mechanisms  in  single  phase, 

dual  phase  and  fully  lamellar  titanium  aluminides  using  finite 
element  analysis. 

The  aim  of  these  numerical  study  efforts  was  to  provide  insight  into  the 
mechanisms  governing  spray  processing  and  mechanical  behavior;  hence  to 
optimize  processing  parameters,  microstructures,  and  mechanical  properties. 


SUMMARY 


The  principal  developments  during  the  three-year  contract  period  may  be 
described  as  follows:  design  and  successful  implementation  of  spray  forming 
experiments  using  y-TiAl  alloys;  microstructural  characterization  and 
mechanical  testing  on  spray  formed  y-TiAl  alloys;  two  dimensional  modeling 
of  momentum  and  thermal  behavior  during  spray  atomization  of  y-TiAl; 
investigation  of  thermal  residual  stresses  in  spray  formed  y-TiAl  using  finite 
element  methods;  and  finite  element  simulation  of  creep  deformation  and 
rupture  of  titanium  aluminides.  These  developments  are  summarized  in  the 
following  paragraphs. 

Gamma  titanium  aluminides  with  three  different  alloy  compositions  (Ti- 
37.4Al-1.5Nb-1.5Cr,  Ti-42Al-2Nb-2Cr,  Ti-47A1,  all  in  atomic  percent)  were 
successfully  spray  formed.  The  spray  formed  y-TiAl  deposits  generally 
exhibited  a  Gaussian  profile.  The  spray  formed  Ti-47at.%Al  alloys  generally 
exhibited  a  fine,  equiaxed  fully  lamellar  a2+y  structure  (an  average  linear 
grain  size  of  60  pm  and  an  average  interlamellar  spacing  about  0.3  pm),  in 
contrast  with  the  coarse  grains  (typically  250  to  1800  pm)  of  y-TiAl  with  fully 
lamellar  microstructures  that  are  produced  from  conventional  cast  or  cast  + 
heat  treatment  processes.  This  is  significant,  since  it  has  been  reported  that  a 
decrease  in  lamellar  grain  size  increases  both  strength  and  ductility  of  y-TiAl 
with  fully  lamellar  microstructures.  It  is  anticipated  that  a  fully  lamellar  y- 
TiAl  with  a  fine  grain  size  will  lead  to  a  good  combination  of  ductility, 
strength,  fracture  toughness,  and  creep  resistance.  Similar  to  Ti-47A1,  the 
grain  refinement  through  spray  forming  processing  was  also  significant  for 
Ti-37.4Al-l.5Nb-l.5Cr  and  Ti-42Al-2Nb-2Cr.  For  Ti-37.4Al-l.5Nb-l.5Cr,  the 
average  linear  size  of  prior  P  phase  grains  had  been  reduced  from  more  than 
1000  pm  to  approximate  146  pm,  while  that  for  Ti-42Al-2Nb-2Cr  was  reduced 
from  425  pm  to  130  pm.  Other  microstructural  details,  however,  were 
dramatically  different  than  that  in  Ti-47A1.  In  contrast  with  the  fully  lamellar 
a2+y  structure  observed  in  Ti-47A1,  no  lamellar  a2+ywere  detected  in  Ti- 
37.4Al-1.5Nb-1.5Cr.  Rather,  Ti-37.4Al-l.5Nb-l.5Cr  consists  of  prior  P  grain 
boundaries,  with  a  plates  precipitated  out  of  P  phase  in  the  grain  along 
several  preferred  directions,  forming  Widmanstatten  structures.  The  a  plates 


ordered  into  0.2  upon  further  cooling.  The  microstructure  of  Ti-42Al-2Nb- 
2Cr,  on  the  other  hand,  is  featured  by  equiaxed  fully  lamellar  a2+y  structure, 
with  the  02+7  lamellae  being  cross-cut  by  a  series  of  nearly  parallel  lines  in 
each  grain. 

The  compressive  creep  behavior  of  spray  formed  y-TiAl  with  a  fine,  equiaxed 
fully  lamellar  (FL)  microstructure  was  studied  in  a  temperature-stress  regime 
of  780  to  850  °C  and  180  to  320  MPa.  An  apparent  stress  exponent  of  4.3  and  an 
activation  energy  of  342  kj/mol  were  observed  in  the  high-temperature  high- 
stress  regime.  Compared  with  the  FL  y-TiAl  which  was  obtained  through 
conventional  casting  +  heat  treatment  processes,  the  spray  formed  y-TiAl 
exhibited  higher  creep  resistance.  The  higher  creep  resistance  observed  in  the 
present  study  was  discussed  in  light  of  the  interstitial  level,  the  chemical 
composition,  the  grain  size,  and  the  interlocking  of  lamellae  at  the  grain 
boundary  which  in  turn  may  be  a  function  of  interlamellar  spacing  and  the 
step  height  of  the  serrated  grain  boundaries.  It  was  suggested  that  the  small 
interlamellar  spacing  and  possibly  larger  step  height  may  contribute  to  the 
higher  creep  resistance  observed  in  the  present  study. 

A  two-dimensional  model  was  formulated  in  the  study  of  the  momentum 
and  thermal  behavior  during  spray  atomization  of  y-TiAl.  The  modeling 
results  indicated  that  the  velocity,  temperature,  cooling  rate,  flight  time  and 
solidification  behavior  of  atomized  droplets  strongly  depends  on  the  initial 
position,  ro,  and  the  diameter  of  the  droplet,  D.  In  addition,  the  velocity 
profiles,  temperature  histories,  cooling  rates,  flight  times  and  solidification 
behavior  of  droplets  are  also  closely  related  to  each  other.  Specifically,  for 
identical  ro,  i)  fine  droplets  are  more  readily  accelerated  or  decelerated  relative 
to  coarse  droplets  along  both  axial  and  radial  directions,  ii)  fine  droplets  cool 
at  a  higher  cooling  rate  than  coarse  ones,  and  iii)  fine  droplets  also  solidify  at  a 
faster  rate  relative  to  coarse  ones.  For  identical  D,  i)  the  maximum  value  of 
the  axial  component  of  droplet  velocity  decreases  with  increasing  ro,  while 
the  maximum  value  of  the  radial  component  of  droplet  velocity  initially 
increases  and  then  decreases  with  increasing  ro,  ii)  flight  time  increases  with 
increasing  ro  and  iii)  the  solid  fraction  at  any  given  axial  distance  increases 
with  increasing  ro  resulting  from  the  longer  flight  time  for  larger  ro. 
Moreover,  the  two  dimensional  droplet  size  distribution  in  the  spray  cone 


was  demonstrate  to  change  from  being  heterogeneous  to  be  almost 
homogeneous  with  increasing  axial  flight  distance.  When  the  axial  distance 
is  smaller  than  a  critical  value  (approximately  0.75  m  for  the  conditions  used 
in  present  study),  coarse  droplets  concentrate  in  the  central  region  of  the  spray 
cone,  whereas  fine  droplets  populate  the  periphery.  When  the  axial  distance 
increases  beyond  this  critical  value,  the  spatial  droplet  size  distribution  in  the 
spray  cone  becomes  independent  of  the  radial  position.  Finally,  the  two 
dimensional  distribution  of  the  fraction  solidified  in  the  spray  cone  was  noted 
to  be  heterogeneous.  The  fraction  solidified  at  any  axial  distance  increases 
with  increasing  radial  distance  from  the  spray  axis. 

Thermo-elastoplastic  finite  element  analysis  was  carried  out  to  investigate  the 
thermal  residual  stresses  that  may  develop  during  the  spray  deposition 
processing  of  y-TiAl.  It  was  noted  that  thermal  residual  stresses  were 
developed  in  the  spray  formed  y-TiAl  during  cooling  from  the  deposition 
temperature  to  ambient  temperature.  The  axial,  hoop,  and  radial  stresses 
were  found  to  be  mostly  compressive  for  the  lower  center  region  of  the 
preform.  The  magnitude  of  these  compressive  stresses  increases  for  the 
region  closer  to  the  bottom  of  the  preform.  For  the  other  region,  the  axial, 
hoop,  and  radial  stresses  become  tensile.  In  addition,  the  residual  stress 
distribution  was  found  to  be  a  strong  function  of  cooling  rate  in  the  deposited 
preform,  with  the  residual  stresses  decrease  as  the  cooling  rate  decreases. 
When  the  substrate  is  insulated,  the  preform  has  the  lowest  residual  stress 
distribution  for  the  cases  studied  herein.  Finally,  the  residual  stress 
distribution  was  also  demonstrated  to  be  a  strong  function  of  deposit 
thickness.  For  the  region  closer  to  the  outer  edge  of  the  preform,  which 
experiences  the  most  severe  thermal  gradient,  has  the  highest  residual 
stresses.  The  von  Mises'  equivalent  stresses  in  this  region  were  found  to  be 
greater  than  the  yield  stress  of  y-TiAl.  Therefore,  it  is  anticipated  that  plastic 
deformation  or  cracks  may  be  present  in  this  region  of  the  preform. 

Creep  deformation  and  cavitation  damage  in  dual  phase  (a2+y)  equiaxed  and 
fully  lamellar  microstructures  were  simulated  and  analyzed  in  the  present 
research  program  using  finite  element  techniques.  Different  models  of  these 
microstructures  were  developed  depending  on  the  different  deformation 
mechanisms  active  in  the  microstructures  such  as  grain  boundary  sliding  and 


phase  boundary  sliding.  Nonlinear  viscous  primary  creep  deformation  was 
modeled  in  each  phase  based  on  published  creep  data.  Overall  strain  rates  are 
compared  to  gain  an  understanding  of  the  relative  influence  each  of  these 
localized  deformation  mechanisms  has  on  the  creep  strength  of  the 
microstructures  considered.  Facet  stress  enhancement  factors  were  also 
determined  for  the  transverse  grain  facets  in  each  model  to  examine  the 
relative  susceptibility  to  creep  damage.  Creep  constrained  grain  boundary 
cavitation  was  also  analyzed  for  these  microstructures.  The  cavitation  in  the 
models  was  based  on  the  modified  equations  of  Rice  and  Needleman.  It  was 
found  that  grain  boundary  sliding  strongly  enhances  the  cavity  growth  in  all 
of  the  models  analyzed.  Comparisons  of  the  results  for  different  lath 
configurations  suggest  that  lath  width  has  little  effect  on  the  constraints  that 
govern  cavity  growth.  A  study  of  effect  of  cavitating  facet  density  on 
cavitation  damage  was  completed  for  the  fully  lamellar  microstructures. 
Approximate  techniques  of  Hutchinson  and  Riedel  and  Eggeler  for  strain  rate 
enhancement  during  creep  cavitation  due  to  increase  in  cavitating  facet 
density  were  compared  using  the  finite  element  results  as  the  reference. 
Interaction  between  cavitating  facets  in  the  fully  lamellar  TiAl  were  been 
studied.  It  was  found  that  when  the  cavitating  facets  are  not  on  neighboring 
grain  boundaries,  a  negative  interaction  effect  takes  place  in  the  presence  of 
grain  boundary  sliding.  The  present  results  suggest  that  the  longer  creep  life 
observed  experimentally  for  the  fully  lamellar  structure  is  primarily  due  to 
inhibited  former  grain  boundary  sliding  in  this  microstructure  compared  to 
relatively  unimpeded  grain  boundary  sliding  in  the  equiaxed  microstructures. 
The  serrated  nature  of  the  former  grain  boundaries  generally  observed  for 
lamellar  TiAl  alloys  is  consistent  with  this  finding. 
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CHAPTER  1. 


BACKGROUND 


The  ordered  y-TiAl  intermetallic  is  an  attractive  candidate  material  for 
applications  that  demand  thermal  stability,  as  a  result  of  its  attractive 
combinations  of  density  (3.76  g/cm^),  creep/stress  rupture  and  fatigue 
resistance,  and  elevated  temperature  strength  [1-10].  Not  surprisingly,  the 
available  literature  reveals  that  both  single  phase  and  two  phase  y-TiAl  alloys 
have  been  studied  extensively  [11-14].  The  physical  and  mechanical 
properties  of  multicomponent  y-TiAl  alloys  have  been  demonstrated  to 
depend  on  chemistry,  phase  composition  and  geometry,  thermomechanical 
processing,  and  the  content  of  interstitial  impurities  (e.g.,  nitrogen,  carbon, 
oxygen  and  boron)  [12-17].  Extensive  structural  utilization  of  y-TiAl, 
however,  has  been  hindered  by  the  poor  low  temperature  ductility  and  less- 
than-optimum  toughness  that  are  characteristic  of  this  material.  In 
particular,  the  solute  segregation  (both  macrosegregation  and 
microsegregation)  that  is  generally  associated  with  cast  y-TiAl  requires 
extensive  post  processing  long  term  heat  treatments,  thereby  leading  to  low 
material  yields  and  high  costs  [2-4,  11,  18-20].  Furthermore,  although 
macrosegregation  may  be  eliminated  subsequently  by  either  hot  working  (e.g., 
extrusion,  forging,  hot  isostatic  pressing),  heat  treatment  (typically  in  the 
alpha  field  for  cast  y-TiAl)  or  both  [11],  the  suppression  of  microsegregation 
requires  an  enhancement  of  solidification  kinetics  during  processing.  To  that 
effect,  it  is  possible  to  avoid  macrosegregation  and  minimize 
microsegregation  in  y-TiAl  by  using  powder  metallurgy  techniques. 
Nonetheless,  powder  metallurgy  technique  inherently  introduces  surface 
contaminants  (as  a  result  of  the  large  surface  areas  that  are  associated  with 
fine  powders),  and  leads  to  higher  manufacturing  costs.  Inspection  of  the 
available  scientific  literature  reveals  that  both  intrinsic  (e.g.,  manipulation  of 
grain  morphology,  alloying  additions  and  novel  processing  methods)  and 
extrinsic  (e.g.,  ductile  phase  toughening)  approaches  are  actively  being 
pursued  in  an  effort  to  improve  the  properties  of  y-TiAl  [11,  13,  14,  17].  It  is 
also  evident,  however,  that  existing  processing  methods  should  be  optimized 


and  novel  processing  methods  should  be  explored,  before  y-TiAl  can  be 
successfully  exploited  as  a  structural  material. 

In  view  of  the  above  findings,  present  research  efforts  are  aimed  at  using  a 
novel  synthesis  approach,  namely  spray  atomization  and  deposition,  to 
eliminate  macrosegregation,  minimize  microsegregation,  and  achieve  grain 
refinement  and  compositional  homogeneity  in  y-TiAl.  In  spray  atomization 
and  deposition,  a  molten  stream  of  metal  is  disintegrated  into  a  fine 
dispersion  of  droplets  using  high  velocity  inert  gas  jets.  The  resulting  semi- 
solidified  droplets  are  directed  towards  a  substrate  where  they  impinge  and 
collect  as  fine  grained,  chemically  homogeneous,  thin  splats.  The  first  high 
velocity  droplets  impact  directly  against  the  substrate  surface,  and  then 
themselves  act  as  a  quenching  substrate  for  their  successors.  Relatively  high 
rates  of  solidification  are  imposed  during  processing  as  a  result  of  the 
thinness  of  the  splats,  the  rapid  heat  extraction  during  flight,  and  an  inert 
atmosphere.  Spray  atomization  and  deposition  inherently  avoids 
microstructural  coarsening  and  extensive  macrosegregation,  which  are 
normally  associated  with  casting  processes.  To  that  effect  it  is  well 
documented  that  spray  atomization  and  deposition  has  been  successfully 
applied  to  improve  the  physical  and  mechanical  characteristics  of:  A1  and  Mg 
alloys  [21-23];  Ni,  Fe  and  NisAl  intermetallics  [24-27];  and  discontinuous 
metal  matrix  composites  (MMCs)  [28-31].  More  recently,  encouraging 
preliminary  results  were  reported  on  spray  atomization  and  deposition 
processing  of  TiAl  and  Ti^Al  [32-35].  J.W.  Sears  [34],  for  example,  notes  that 
initial  spray  deposition  experiments  led  to  low  interstitial  levels,  low 
porosity,  and  fine  grained  structures.  The  latter  observation,  namely,  a  fine 
grained  structure,  should  lead  to  improvements  in  toughness.  In  other 
studies  Sears  et  al.  [35]  were  able  to  produce  a  totally  dense  material  using  Hot 
Isostatic  Pressing  (HIPing)  following  spray  deposition.  J.M.  Young  et  al.  [33] 
also  report  that  the  pores  that  were  present  in  the  as-sprayed  Ti-52A1  were 
removed  by  HIPing.  None  of  these  studies,  however,  reported  any  processing 
details.  Moreover,  evaluation  of  the  microstructural 

characteristics/variations  and  mechanical  properties  are  limited  in  these 
studies,  not  to  mention  correlation  of  processing  parameter  with 
microstructure/mechanical  properties  of  spray  formed  y-TiAl. 


CHAPTER  2. 


SPRAY  FORMING  OF  y-TiAl  ALLOYS 

2.1.  Objective 

The  efforts  in  this  aspect  include  the  design,  and  the  implementation,  of  spray 
forming  experiments  using  y-TiAl  alloys;  and  characterization  of  the  spray 
formed  alloys. 
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2.2.  Experimental 

The  induction  skull 
melting/spray  atomization  and 
deposition  facility  utilized  to 
spray  process  gamma  TiAl 
alloys  is  schematically  shown 
in  Figure  2-1.  The 
experimental  procedure  in 
spray  forming  may  be  described 
as  follows.  The  skull  crucible 
was  charged  with  as-received  y- 
TiAl  alloys.  The  cooling 
system  in  the  skull  crucible 
assured  the  formation  of  a  y- 
TiAl  skull  on  the  inside  wall  of 

the  crucible  when  the  molten  Schematic  representation  of  the 

induction  skull  melting/spray  atomization  and 
Ti-Al  alloys  contact  the  cold  deposition  facility. 

crucible  wall.  The  formation  of 

the  y-TiAl  skull  prevented  the  reaction  between  the  chemically  reactive 
molten  Ti-Al  alloys  and  the  crucible.  Once  the  central  part  of  the  charge  was 
completely  melted,  the  skull  crucible  was  tilted  down,  and  Ti-Al  melt  was 
poured  into  a  synchronously  heated  tundish.  The  Ti-Al  melt  was  then 
delivered  through  a  nozzle  to  an  18-jets  atomizer,  where  it  was  atomized  by 
high  pressure  Argon  gas  (dynamic  atomization  pressure  ranges  from  1  to  2 
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MPa).  The  atomized  micro-sized  droplets  were  subsequently  deposited  on  a 

water  cooled  stationary  planar  Cu  substrate.  Powders  not  incorporated 

into  the  deposit,  which  are  generally 

referred  to  as  oversprayed  powders, 

were  carried  by  the  gas  flow  into  the 

cyclone  and  collected  in  the 

container  attached  to  the  cyclone. 

To  prevent  extensive  oxidation  of 

the  materials  during  processing,  all 

experiments  were  conducted  in  an 

environmental  chamber,  which  was 

„  ,  ,  Figure  2-2.  Gaussian  profile  of  spray 

generally  evacuated  to  a  pressure  of  formed  Ti-47at.%Al. 

133  Pa  (1  torr)  or  lower,  and  back 

filled  with  argon  gas  to  a  pressure  of  1x10'^  Pa  prior  to  melting  and 
atomization. 


Gamma  titanium  aluminides  with  three  different  alloy  compositions  were 
successfully  spray  formed,  as  summarized  in  Table  2-1.  The  spray  formed  y- 
TiAl  deposits  generally  exhibited  a  Gaussian  profile,  as  shown  in  Figure  2-2 
for  a  spray  formed  Ti-47A1.  The  thickness  of  the  central  region  ranges  from  10 
to  15  mm  depending  on  processing  parameters,  while  the  diameter  of  the 
bottom  ranges  from  10  to  15  cm. 


Table  2-1.  Spray  formed  y-TiAl  alloys 


Nominal 

composition  (at.%) 

Supplier 

Ti-37.4Al-l.5Nb- 

1.5Cr 

Supplied  by  Dr.  Don  Larsen,  Howmet  Corporation, 
Operhall  Research  Center,  Whitehall,  MI 

Ti-42Al-2Nb-2Cr 

Supplied  by  Dr.  Shyh-Chin  Huang,  General  Electric 
Research  and  Development  Center,  Schenectady,  NY 

Ti-47A1 

Purchased  from  The  Duriron  Company,  Inc., 
Titanium  Casting  Operations,  Dayton,  OH 

2.3.  Microstructures 

The  spray  formed  Ti-47at.%Al  alloys  generally  exhibited  a  fine,  equiaxed  fully 
lamellar  a2+y  structure.  Figure  2-3  shows  the  microstructure  of  one  of  the 
spray  formed  Ti-47at.%Al,  with  an  average  linear  grain  size  of  60  |im  and  an 


average  interlamellar  spacing  about  0.3  pm,  in  contrast  with  the  coarse  grains 
(typically  250  to  1800  pm  [5,  36])  of  y-TiAl  with  fully  lamellar  microstructures 
that  are  produced  from  conventional  cast  or  cast  +  heat  treatment  processes. 
This  is  important,  since  it  has  been  reported  that  a  decrease  in  lamellar  grain 
size  increases  both  strength  and  ductility  of  y-TiAl  with  fully  lamellar 
microstructures  [5].  It  is  anticipated  that  a  fully  lamellar  y-TiAl  with  a  fine 
grain  size  will  lead  to  a  good  combination  of  ductility,  strength,  fracture 


Figure  2-3.  Optical  microstructure  of  spray  formed  Ti-47at.%Al. 


The  spray  formed  Ti-37.4Al-l.5Nb-l.5Cr  exhibits  an  equiaxed  morphology, 
with  an  average  linear  grain  size  of  approximately  146  pm,  as  shown  in 
Figure  2-4.  The  acicular  a  phase  precipitates  from  the  elevated  temperature  P 
phase,  forming  a  Widmanstatten  structure.  The  grain  refinement  of  Ti- 
37.4Al-1.5Nb-1.5Cr  through  spray  forming  processes  is  significant,  since  the 
grain  size  of  the  as-received  material  is  larger  than  1000  pm.  In  comparison, 
spray  formed  Ti-42Al-2Nb-2Cr  exhibits  a  fully  lamellar  a2+y  microstructure, 
with  an  average  grain  size  of  approximately  130  pm  (Figure  2-5),  refined  from 
an  average  grain  size  of  425  pm  for  the  as-received  materials.  It  is  worth 
noting  that,  in  each  grain,  there  are  several  nearly  parallel  lines  across  the 
a2+y  lamellae,  as  marked  in  Figure  2-5  in  one  of  the  grains.  Based  on  the 
results  of  heat  treatment  on  as-received  and  spray  formed  Ti-37.4Al-l.5Nb- 


1.5Cr  (not  Ti-42Al-2Nb-2Cr),  it  is  believed  that  these  lines  are  the  traces  of 
acicular  a  phase  grain  boundaries  (aPGB),  i.e.,  lamellar  a2+Ywere 


Figure  2-4.  Optical  microstructure  of  spray  formed  Ti-37.4Al-l.5Nb-l.5Cr. 


Figure  2-5.  Optical  microstmcture  of  spray  formed  Ti-42Al-2Nb- 
2Cr.  aPGB  represents  prior  acicular  a  phase  grain  boundaries. 
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2.4.  Summary 

Gamma  titanium  aluminides  with  three  different  alloy  compositions  were 
successfully  spray  formed.  Significant  structure  refinement  through  spray 
atomization  and  deposition  process  was  observed  for  all  these  three  alloys. 


CHAPTER  3. 


COMPRESSIVE  CREEP  BEHAVIOR  OF  SPRAY  FORMED 
GAMMA  TITANIUM  ALUMINIDE 


3.1.  Objective 

The  purpose  of  this  effort  is  to  investigate  the  compressive  creep  behavior  of 
spray  formed  y-TiAl,  and  compare  it  with  that  of  conventionally  processed  y- 
TiAl. 


3.2.  Experimental 

Rectangular  parallelepiped  compression  specimens  with  a  3  mm  gauge  length 
and  a  1.5  mm  x  1.5  mm  cross  section  were  prepared  from  the  spray  formed  Ti- 
47A1  alloy.  Uniaxial  compression  creep  tests  were  conducted  in  air  over  a 
temperature  range  of  780-850  °C  (0.58-0.62  Tm,  where  T^  is  the  absolute 
melting  temperature  of  y-TiAl)  under  true  stresses  ranging  from  180  to  320 
MPa.  All  tests  were  performed  at  a  single  constant  true  stress  and 
temperature.  The  microstructure  prior  and  following  creep  deformation  was 
examined  using  optical  microscopy.  The  etchant  used  was  Kroll's  reagent 
with  the  following  composition:  1  part  HF,  2  parts  HNO3  and  50  parts  H2O. 


Figure  3-1.  Creep  curve  of  spray  formed  y-TiAl  at  850  °C  with  a  stress  of  300  MPa: 

(a)  strain  as  a  function  of  time,  and  (b)  creep  rate  as  a  function  of  strain. 

3.3.  Results 

Figure  3-1  shows  a  typical  creep  curve  obtained  for  spray  formed  y-TiAl.  The 
data  is  for  a  temperature  of  850  °C  and  a  true  stress  of  300  MPa.  A  primary 


creep  stage  (where  the  creep  rate  decreases  with  strain)  was  observed  before  a 
minimum  creep  rate  of  1.04  x  lO'^  S'^  was  attained.  This  was  followed  by  a 
region  in  which  the  creep  rate  gradually  increases  with  strain.  The 
experimentally  observed  minimum  creep  rates  of  all  tests  in  the  present  study 
are  summarized  in  Table  3-1. 


Table  3-1.  Summary  of  compressive 
creep  testing  of  spray  formed  y-TiAl 


Temperature  (°C) 

True  stress  (MPa) 

Minimum  creep  rate  (s‘0 

300 

1.04x10-6 

260 

5.76  X  10-7 

220 

3.25  X  10-7 

180 

8.57  X  10-8 

825 

300 

3.69  X  10-7 

825 

280 

3.1  X  10-7 

320 

1.71  X  10-7 

800 

300 

1.66x10-7 

280 

8.43  X  10-8 

800 

250 

5.73  X  10-8 

300 

8.97  X  10-8 

In  general,  the  creep  deformation  of 
phenomenological  equation  such  as: 

e  min  =  A  exp(-Q/RT) 
where  Emin  is  the  minimum  strain 
rate,  A  is  a  pre-exponential  constant, 
a  is  the  applied  stress,  n  is  the  stress 
exponent,  Q  is  the  activation  energy 
for  creep,  R  is  the  universal  gas 
constant,  and  T  is  the  absolute 
temperature.  Figure  3-2  shows  the 
variation  of  minimum  creep  rate  as  a 
function  of  stress  at  temperatures  of 
800  and  850  °C.  The  slope  of  the 
curves  in  Figure  3-2  yields  the  stress 
exponent,  n,  according  to  Eqn.  (3-1). 
From  Figure  3-2,  the  values  of  the 
stress  exponent,  n,  were  identified  to 


a  material  may  be  described  by  a 


(3-1) 


100  1000 
Stress  (MPa) 


Figure  3-2.  Minimum  creep  rate  as  a 
function  of  stress  at  800  °C  and  850  °C. 


be  4.9  and  4.8  for  800  °C  and  850  °C,  respectively.  Figure  3-3  shows  the 
variation  of  minimum  creep  rate  as  a  function  of  temperature  with  a  stress  of 
300  MPa.  An  activation  energy  for  creep  of  342  kj/mol  was  determined  from 
this  plot.  Figure  3-4  shows  the  temperature  compensated  minimum  creep 
rate  as  a  function  of  stress  applied,  using  an  activation  energy  for  creep  of  342 
kJ/mol,  for  data  shown  in  Figure  3-2  plus  additional  data  at  780  and  825  °C. 
All  creep  data  obtained  under  different  temperatures  and  stresses  cluster 
about  a  straight  line,  as  predicted  from  Eqn.  (3-1),  with  a  value  of  the  slope 
equal  to  4.3.  Curve  fitting  to  the  data  shown  yields  the  following  equation  to 
describe  the  creep  behavior  of  spray  formed  y-TiAl: 

e  min  =  0.21  04.3  exp(-342kJmol-i  /RT)  (3-2) 


Figure  3-3.  Minimum  creep  rate  as  a  Figure  3-4.  Temperature  compensated 
function  of  temperature  at  300  MPa.  minimum  creep  rate  as  a  function  of  stress. 


3.4.  Discussion 

The  observed  increase  in  creep  rate  with  strain  in  Figure  3-1  at  strains  greater 
than  4.5%  may  imply  that  a  microstructural  change  occurs  in  the  specimen 
during  deformation.  To  that  effect,  it  has  been  suggested  that  a 
spheroidization  of  the  lamellar  microstructure  may  occur  under  thermo¬ 
mechanical  exposure  [37-39].  However,  as  shown  in  Figure  3-5,  optical 
microstructural  observation  of  the  y-TiAl  following  creep  testing  at  850  °C 
and  300  MPa  revealed  no  obvious  spheroidization.  Other  factors  which  can 
lead  to  the  increase  in  creep  rate  with  strain  include:  (1)  void  formation  or 
loss  of  internal  section,  (2)  environmental  degradation  (i.e.,  oxidation),  (3) 


necking  or  loss  of  external  section  which  is  normally  associated  with  tensile 
testing,  (4)  microstructure  degradation  (i.e.,  particle  coarsening  or  grain 
growth),  and  (5)  inverse  creep  which  is  associated  with  the  glide  of 
dislocations  [40-42].  Comparison  between  optical  microstructures  of  as  spray 
formed  and  crept  specimens  reveals  no  obvious  grain  growth  during  the 
creep  testing.  Moreover,  as  will  be  shown  below,  it  appears  that  the 
governing  mechanism  for  the  creep  of  spray  formed  y-TiAl  is  dislocation 
climb.  Therefore,  for  the  compressive  creep  behavior  of  spray  formed  y-TiAl, 
it  is  unlikely  that  the  observed  increase  in  creep  rate  with  strain  results  from 


either  of  the  last  three  factors  mentioned  above.  Conclusive  determination  of 


Figure  3-5.  Microstructure  of  y-TiAl  following  creep  testing,  8=12.11%,  at 
850  °C  with  a  stress  of  300  MPa.  The  arrow  indicates  the  distorted  lamellae. 


the  effect  of  environmental  degradation  requires  a  comparison  of  creep 
behavior  in  vacuum  and  in  air.  However,  since  little  oxidation  has  been 
observed  on  the  surface  of  the  crept  specimens  in  the  present  study,  it  is  likely 
that  the  increase  of  creep  rate  with  strain,  after  the  minimum  creep  rate  has 
been  obtained,  results  from  the  formation  of  voids.  The  voids  in  specimen 
crept  at  850  °C  with  a  stress  of  300  MPa  until  8=12.11%  is  shown  Figure  3-6a. 
Since  there  are  also  voids  in  as-spray  formed  materials,  as  shown  in  Figure  3- 
6b,  more  reliable  evidence  on  formation  of  voids  during  creep  testing  should 
be  the  general  trend  of  density  decrease  for  each  specimen  after  creep 


deformation.  This  necessitates  the  measurement  of  density  changes  for  all  of 
the  specimens  prior  to  and  after  the  creep  deformation,  which  has  not  been 
carried  out  in  the  present  study.  It  is  worth  noting  that  in  a  recent  study 
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Figure  3-6.  Voids  observed  under  optical  microscopy  on  metallographically 
polished/unetched  y-TiAl  specimens:  (a)  crept,  and  (b)  spray  formed. 


conducted  by  Es-Souni  et  al.  [8]  on  creep  behavior  of  fully  lamellar  y-TiAl 
obtained  through  conventional  cast  +  heat  treatment,  TEM  observation 


revealed  no  dynamic  recrystallization  of  single  phase  y  regions  or 
spheroidization  of  02  laths,  which  was  detected  to  be  the  fundamental  reasons 
for  the  observed  minimum  creep  rate  in  a  y-TiAl  with  lamellar  +  y 
microstructure  [43].  Accordingly,  early  formation  of  voids  and  cavities  at 
colony  boundaries  due  to  stress  concentration,  which  may  result  from  the 
strong  dependence  of  yield  stress  on  the  orientation  of  the  lamellae  with 
respect  to  the  tensile  stress,  was  suggested  to  be  responsible  for  the  absence  of 
steady  state  in  creep  behavior  of  fully  lamellar  y-TiAl  [8]. 

A  stress  exponent  of  4.3  in  the  present  study  indicates  that  the  creep  behavior 
of  the  spray  formed  y-TiAl  is  likely  controlled  by  a  dislocation  climb  process 
[44-46].  In  this  case,  the  activation  energy  for  creep  should  be  equal  to  the 
activation  energy  for  lattice  diffusion  of  the  slowest  species  (Ti  or  Al)  in  y- 
TiAl.  Unfortunately,  little  data  is  available  on  self-diffusion  of  Ti  or  Al  in 
two  phase  y-TiAl.  It  is  worth  noting,  however,  that  the  activation  energy  for 
self-diffusion  of  Ti  in  single  phase  y-TiAl  (Ti-54A1)  has  been  measured  to  be 
291  kj/mol  [47].  This  is  close  to  the  activation  energy  for  creep,  342  kj/mol,  as 
observed  in  the  present  study. 

The  creep  behavior  of  y-TiAl  based  alloys  with  different  microstructures 
(fully  lamellar,  duplex,  nearly  lamellar,  and  single  phase  gamma)  has  been 
investigated  extensively,  and  turns  out  to  be  complex.  Even  for  fully  lamellar 
y-TiAl  alone,  inspection  of  the  available  literature  indicates  that  there  are 
some  discrepancies  in  the  reported  data.  While  creep  curves  similar  to  that 
shown  in  Figure  3-1  have  been  observed  by  several  other  researchers  [8,  38, 
48],  a  well  defined  steady  state  was  also  reported  by  Worth  et  al.  [49]  and  Es- 
Souni  et  al.  [50]  for  fully  lamellar  y-TiAl.  This  may  result  from  the 
differences  in  alloy  composition  or  microstructural  variables,  such  as  grain 
size  and  interlamellar  spacing,  associated  with  each  specimen  under 
investigation.  The  magnitudes  of  stress  exponent,  n,  and  activation  energy, 
Q,  obtained  in  the  present  study  are  in  agreement  with  that  reported  by  Hayes 
et  al.  (n=4.6  and  Q=335  kJ/mol)  [48],  Bartholomeusz  et  al.  (n=5.4  and  Q=407 
kJ/mol)  [38],  and  Worth  et  al.  (n=5.2,  Q=348  kJ/mol  for  high  stress  level)  [49]. 
As  mentioned  before,  a  stress  exponent  ranging  from  4  to  5  and  an  activation 
energy  for  creep  close  to  that  for  self-diffusion  suggest  a  dislocation  climb 
mechanism  in  creep  deformation.  Es-Souni  et  al.  [8,  50]  have  reported  a 


higher  stress  exponent  (n=  7.6)  with  microstructural  evidence  for  dislocation 
climb  (i.e.  subgrain)  for  fully  transformed  lamellar  y-TiAl.  No  activation 
energy  for  creep  was  determined  in  their  studies.  Thus,  the  present  results 
are  in  agreement  with  previous  results,  including  those  of  Es-Souni  et  al  [8], 
on  that  the  creep  rate  of  fully  lamellar  y-TiAl  is  controlled  by  dislocation 
climb.  A  stress  exponent  of  7.6,  which  was  reported  by  Es-Souni  et  al  [8,  50], 
may  be  rationalized  in  terms  of  two  possible  mechanisms:  1)  dislocation  climb 
under  constant  structure  [51]  or  2)  the  presence  of  a  threshold  stress  [10].  The 
high  stress  exponent  may  result  from  the  presence  of  fine  intermetallic 
dispersoids  or  oxides /nitrides /carbides  in  the  specimen.  The  fine  particles 
can  pin  the  subgrain  boundaries,  stabilizing  the  subgrain  size  (constant 
structure),  or  if  the  spacing  between  particles  is  smaller  than  the  subgrain  size, 
they  can  act  as  effective  barriers  to  dislocation  motion,  giving  rise  to  a 
threshold  stress.  It  should  be  noted  that  for  both  creep  under  constant 
structure  and  creep  in  the  presence  of  a  threshold  stress,  the  creep  rate  can  still 
be  controlled  by  dislocation  climb. 


It  would  be  of  interest  to  compare 
the  minimum  creep  rates 
observed  in  the  present  study  to 
those  reported  by  other  authors. 
To  that  effect,  temperature 
compensated  creep  rate,  using  an 
activation  energy  for  creep 
averaged  on  reported  values  by 
Bartholomeusz  et  al.  [38],  Hayes  et 
al.  [48],  Worth  et  al.  [49],  and  the 
present  study,  was  plotted  versus 
applied  stress,  as  shown  in  Figure 
3-7.  A  stress  exponent  of  4.8, 
averaged  based  on  those  authors. 
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Figure  3-7.  Comparison  of  temperature 
compensated  minimum  creep  rates  as  a  function 
was  used  to  curve  fit  the  of  stress. 

experimental  results.  This  yields 
the  following  equations: 

=  0.251  o4-8exp(-358kJmol‘i/RT) 

£  min  =  0.2  ®  exp(-358kJmol'i/RT) 


(3-3) 

(3-4) 


for  Worth  et  al.  [49] 
for  Hayes  et  al.  [48] 


8min  =  0-2  cj‘^'^exp(-358kJmoH/RT)  (3-5)  for  Bartholomeusz 

et  al.  [38] 

and  £  min  =  0.058  exp(-358kJinol-i/RT)  (3-6)  the  present  study 


Evidently  from  Figure  3-7  and  the  above  equations,  the  minimum  creep  rate 
observed  in  the  present  study  is  lower  than  those  reported  by  other  authors 
for  a  stress-temperature  range  where  dislocation  climb  is  the  dominant 
mechanism.  For  FL  microstructures,  the  variables  which  may  influence  the 
creep  resistance  include:  interstitial  level,  chemical  composition,  grain  size 
and  interlamellar  spacing  [5,  37].  These  variables  for  the  materials  used  by  the 
different  authors  are  provided  in  Table  3-2,  with  their  possible  effect  discussed 
individually  below. 


Table  3-2.  Linear  intercept  grain  size,  interlamellar 


Materials 

Grain  size 
(pm) 

Interlamella 
r  spacing 
(pm) 

C 

(at.  %) 

N 

(at.  %) 

O 

(at.  %) 

Source 

Ti-48Al-2Nb- 

2Mn 

150-200 

— 

0.035 

0.018 

■ 

[48] 

Ti-49Al-lV 

850 

0.54 

0.07 

0.01 

0.12 

[49] 

Ti-44A1 

70 

5.1 

— 

— 

— 

[38] 

Ti-47A1 

60 

0.3 

0.25 

0.08 

0.35 

Present 

study 

(A).  Interstitial  level 

It  is  evident  from  Table  3-2  that  the  interstitial  levels  of  carbon  and  oxygen  in 
the  material  used  in  the  present  study  are  higher  relative  to  those 
investigated  by  other  authors.  This  may  lead  to  the  incorrect  conclusion  that 
the  higher  creep  resistance  observed  in  the  present  study  resulted  from  the 
higher  interstitial  content  in  the  material  investigated.  When  the  creep 
deformation  is  controlled  by  dislocation  glide,  the  above  conclusion  may 
sustain,  since  the  movement  of  dislocation  along  the  slip  plane  may  be 
dragged  by  the  solute  atmosphere,  leading  to  decreased  creep  rate  (i.e.,  higher 
creep  resistance)  [52].  However,  as  discussed  earlier,  the  creep  behavior  as 
reported  by  Worth  et  al.  [49],  Hayes  et  al.  [48],  Bartholomeusz  et  al.  [38],  and  the 


present  study  was  controlled  by  dislocation  climb,  not  dislocation  glide.  In 
this  case,  two  solute  strengthening  mechanisms  are  possible:  i).  interaction  of 
solute  with  stacking  faults,  and  ii).  pinning  of  the  subgrain  boundary  by  the 
precipitate  /  dispersoid. 

The  interaction  of  solute  with  stacking  faults  may  decrease  the  dislocation 
mobility  by  increasing  the  stress  for  recombination  of  partial  dislocation 
(hence  decreasing  the  ability  to  activate  cross  slip  or  climb)  [52].  Nevertheless, 
it  was  found  that  the  creep  deformation  structure  of  the  fully  lamellar,  at  a 
temperature-stress  regime  similar  to  the  one  investigated  here,  was 
dominated  by  the  presence  of  1/2<110]  dislocations  [8,  52],  of  which  partial 
formation  is  inhibited  since  it  involves  very  high  stacking  fault  energy  [53]. 
When  partial  formation  is  inhibited,  solute  interaction  with  stacking  faults 
and  the  corresponding  strengthening  effect  vanish  [54].  Therefore,  it  is 
unlikely  that  the  interstitial  will  increase  the  creep  resistance  through 
solute/stacking-fault  interaction. 

When  the  subgrain  boundary  is  pinned  by  the  precipitate/dispersoid,  the 
dislocation  mobility  may  be  reduced  because  of  a  change  from  climb  to  shear 
for  dislocations  to  bypass  the  particles.  This  strengthening  mechanism, 
however,  introduces  a  threshold  stress  and  gives  rise  to  an  experimentally 
determined  stress  exponent  of  8  or  higher  [10,  52,  55,  56].  Indeed,  in  a  study 
conducted  by  Worth  et  al.  [52]  regarding  the  influence  of  carbon  on  creep 
deformation  of  y-TiAl  alloys,  it  was  found  that  the  creep  resistance  of  fully 
lamellar  y-TiAl  was  increased  by  a  factor  of  8  through  increasing  the  carbon 
content  from  0.07  to  0.3  at.%.  Meanwhile,  the  stress  exponent  in  their  study 
[52]  also  increased  from  5.2  for  low  carbon  alloy  to  8.4  for  high  carbon, 
suggesting  that  the  precipitation  of  carbide  in  the  matrix  strengthened  TiAl 
through  a  threshold  stress  mechanism.  In  the  present  study,  the  dominant 
mechanism  for  creep  is  dislocation  climb,  with  the  stress  exponent  in  the 
range  from  4  to  5.  This  excludes  the  influence  of  impurities  on  the  higher 
creep  resistance  observed  in  the  present  study  through  precipitate/ dispersoid 
strengthening  mechanism. 


Therefore,  it  is  proposed  that  the  higher  creep  resistance  observed  in  the 
present  study  is  not  due  to  the  relatively  higher  interstitial  level  associated 
with  the  spray  formed  y-TiAl. 


(B).  Chemical  composition 

The  effect  of  chemical  composition  on  creep  originates  from  its  possible  effect 
on  the  diffusivity.  For  a  NiAl  intermetallic  compound,  it  was  found  that 
when  the  chemical  composition  deviates  from  the  stoichiometric  point,  the 
diffusivity  of  Ni  in  NiAl  increases  [57].  This  phenomenon  was  employed  to 
rationalize  the  observed  increase  of  creep  resistance  with  decreasing  A1 
content  in  single  phase  y-TiAl  [58,  59].  In  the  study  of  single  phase  y-TiAl 
[58,  59],  the  A1  content  is  above  50%  for  all  materials.  Decreasing  A1  content 
corresponds  to  less  deviation  from  the  stoichiometric  composition;  hence 
slower  self-diffusion  and  higher  creep  resistance.  In  the  present  study,  the 
deviation  from  the  stoichiometric  composition  increases  from  Ti-49Al-lV, 
Ti-48-2Nb-2Mn,  Ti-47A1,  to  Ti-44A1.  Accordingly,  it  is  anticipated  that  the 
self-diffusion  of  Ti  or  A1  in  these  materials  should  increase  from  Ti-49Al-lV, 
Ti-48-2Nb-2Mn,  Ti-47A1,  to  Ti-44A1. 

Because  dislocation  climb  is 
essentially  controlled  by  the  lattice 
diffusion,  the  minimum  creep  rate 
is  expected  to  increase  from  Ti-49A1- 
IV,  Ti-48-2Nb-2Mn,  Ti-47A1,  to  Ti- 
44A1.  This  is  completely  different 
from  the  trend  evident  in  Table  3-2, 
suggesting  that  there  are  other 
mechanisms  affecting  the  creep 
resistance. 


(C).  Grain  size 


GL  (urn*) 

Figure  3-8.  Relationship  between  minimum 
From  Figure  3-7,  Equations  3-6,  and  creep  rate,  interlamellar  spacing  (L),  and 

Table  3-2,  it  appears  that  the  higher 

creep  resistance  observed  in  the 

present  study  may  be  associated  with  a  decrease  in  grain  size.  Figure  3-8  plots 
the  pre-exponential  constants  in  Equations  (3),  (5)  and  (6),  which  represent  the 
relative  magnitude  of  creep  rate  for  different  materials  investigated  under  the 


same  temperature  and  stress,  as  a  function  of  the  product  of  grain  size  and 
interlamellar  spacing,  GL.  The  results  in  Figure  3-8  indicate  that  the  change 
of  creep  resistance  may  be  well  correlated  with  the  change  of  GL:  when  GL 
increases,  the  minimum 
creep  rate  increases. 

However,  the  validity  of  this 
type  of  dependence  requires  a 
satisfactory  elucidation  of  its 
fundamental  mechanisms. 

This  is  somehow  difficult  if 
one  notices  that  the  creep 
deformation  is  controlled  by  a 
dislocation  climb  mechanism. 

When  the  creep  deformation 
is  controlled  by  a  dislocation 
climb  mechanism,  it  is 
generally  considered  that 

creep  behavior  is  relatively  insensitive  to  the  change  of  grain  size  [37,  44,  46]. 


Figure  3-9.  TEM  micrograph  showing  the 
interlocking  in  spray  formed  y-TiAl. 


(D).  Interlamellar  spacing 

The  interlamellar  spacing  in  fully  lamellar  microstructure  was  widely 
considered  to  have  major  effect  on  the  creep  behavior  of  y-TiAl.  However,  its 
exact  origin  remains  to  be  explored,  especially  when  the  creep  deformation  is 
controlled  by  dislocation  climb.  It  was  suggested  that  the  increased  creep 
resistance  for  a  lamellar  microstructure  originates  from  the  reinforcing  effect 
of  a2  laths  [12,  52].  Nevertheless,  when  the  deformation  is  controlled  by 
dislocation  climb,  it  is  very  difficult  to  imagine  how  the  a2  laths  act  as 
reinforcement.  More  studies  are  necessary  to  elucidate  its  fundamental 
mechanisms.  Perhaps  it  is  similar  to  the  precipitate/dispersoid 
strengthening.  In  this  case,  a  stress  exponent  of  8  or  higher  should  be 
observed,  which  is  not  true  in  the  present  study. 

The  effect  of  interlamellar  spacing  on  the  creep  behavior  may  be  associated 
with  the  interlocking  of  lamellae  in  neighboring  grains  at  the  grain  boundary. 
The  interlocking  of  lamellae  for  y-TiAl  with  fully  lamellar  microstructure 
was  reported  previously  in  the  literature  by  several  authors  [38,  48,  60].  It  is 


also  observed  in  the  material  used  in  the  present  study,  as  shown  in  Figure  3- 
9.  Presence  of  the  interlocking  may  make  the  grain  boundary  slide  more 
difficult  relative  to  the  case  where  interlocking  is  absent,  as  observed  in 
conventional  materials.  It  has  been  observed  by  Bartholomeusz  et  al.  [38], 
Hayes  et  al.  [48],  as  well  as  being  evident  from  Figure  3-5,  that  the  lamellae  at 
the  grain  boundary  region  are  distorted  after  creep.  This  suggests  that  the 
interlocking  of  lamellae  in  neighboring  grains  at  the  grain  boundary  inhibited 
the  movement  of  grains  relative  to  each  other. 

With  free  grain  boundary  slide,  the  material  may  be  divided  into  soft  and 
hard  bands  parallel  to  the  direction  of  applied  stress  (see  Figure  3-10).  The 
effective  stress  applied  on  the  hard  region  increases  by  a  factor  of  ^  (1<^<2  [61, 
62]),  i.e.,  Ge=^o,  where  Oe  arid  a  are  the  effective  and  applied  stresses, 
respectively.  Without  grain  boundary  slide,  the  polycrystalline  material  in 
Figure  3-10  is  a  rigid  body,  and  the  effective  stress,  Og,  equals  to  the  applied 
stress,  a.  If  the  applied  stress,  a,  in  Eqn.  (3-1),  is  substituted  by  the  effective 
stress  Og,  it  is  evident  that  under  the  same  applied  stress,  the  observed  creep 
rate  would  be  higher  for  the  case  with  grain  boundary  slide  relative  to  that 
without  grain  boundary  slide.  This  concept  was  illustrated  rigorously  in 
reference  [61]  (also  see  Chapter  6),  where  creep  deformation  of  y-TiAl  alloys  at 
two  extreme  cases,  i.e.,  with  and  without  grain  boundary  slide,  were 
investigated.  A  natural  extension  of  this  concept  is  to  incorporate  the 
situations  between  these  two  extreme  cases:  the  effect  of 
smoothness /roughness  of  the  grain  boundary  on  the  effective  stress  on  the 
hard  bands  and,  in  turn,  on  the  creep  resistance  of  the  material.  It  is  not 
difficult  to  envision  that  the  effective  stress  on  the  hard  band  would  be  less 
for  materials  with  rough  grain  boundaries  than  that  with  smooth  grain 
boundaries;  and  the  creep  resistance  of  the  former  would  be  higher  relative  to 
that  of  the  latter.  The  smoothness/roughness  of  the  grain  boundary  may  be  a 
function  of  several  parameters  describing  the  interlocking  of  the  lamellae  of 
neighboring  grains.  First,  it  is  a  function  of  the  interlocking  nodes. 
Assuming  the  interlamellar  spacing  to  be  L,  and  the  segment  of  each  hexagon 
in  Figure  3-10  to  be  S,  total  number  of  nodes  of  interlocking  lamellae  in  each 
interlocking  boundaries  may  be  readily  calculated  to  be  (V3/2)(S/L).  Secondly, 
it  is  intuitively  more  reasonable  to  describe  the  smoothness /roughness  of  the 
grain  boundary  by  the  node  density,  rather  than  the  absolute  number.  The 


node  density  is  given  by  the  total  number  of  nodes  per  unit  length  of  the 
grain  boundary,  which  equals  to  V3/2L.  Finally,  the  roughness  is  also  a 
function  of  the  step  height,  H:  the  larger  the  step  height,  the  more  rough  the 
grain  boundary  is.  Therefore,  the  roughness  of  the  grain  boundary  is 
proportional  to  the  step  height  and  inversely  proportional  to  the 
interlamellar  spacing.  Accordingly,  when  the  step  height  remains  to  be  the 
same,  the  creep  resistance  will  increase  with  decreasing  interlamellar  spacing. 


Figure  3-10.  Schematic  representation  showing:  (a)  the  hard/soft  band 
in  deformation  of  polycrystalline  materials,  and  (b)  the  interlocking 
lamellae  between  neighboring  grain  boundaries.  S,  H,  and  L  are  the 
length  of  the  hexagon  segment,  the  step  height,  and  the  interlamellar 
spacing,  respectively,  s  is  the  applied  stress. 


This  is  consistent  with  the  higher  creep  resistance  observed  in  the  present 
study:  smallest  interlamellar  spacing  corresponds  to  lowest  minimum  creep 

•  rate.  However,  the  minimum  creep  rate  may  not  be  satisfactorily  described  by 
the  interlamellar  spacing  effect  alone  as  evident  from  Figure  3-11,  which 
shows  the  minimum  creep  rate  as  a  function  of  interlamellar  spacing.  This 
suggests  that  the  step  heights  for  the  three  different  materials  in  Figure  3-11 

•  are  different.  Its  effect  has  to  be  taken  into  account.  It  is  worth  noting  that  the 
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chemical  composition  also 
affects  the  creep  resistance, 
as  discussed  earlier. 

Nonetheless,  incorporation 
of  the  chemical 
composition  effect  into  the 
interlamellar  spacing  effect 
does  not  help  to  explain  the 
discrepancy  between  the 

observed  results  and  the  r.  w  u-  u  * 

Figure  3-11.  Relationship  between  minimum  creep 

prediction.  For  example,  rate  and  the  interlamellar  spacing  (L). 
the  creep  resistance 

observed  for  Ti-49Al-lV  was  lower  than  that 
for  Ti-44A1.  According  to  the  interlamellar 
spacing  effect,  the  creep  resistance  for  Ti-49A1- 
IV  should  be  higher  than  that  for  Ti-44A1. 

The  chemical  composition  effect  also  suggests 
that  the  creep  resistance  for  Ti-49Al-lV 
should  be  higher  than  that  for  Ti-44A1. 

Figure  3-12.  Schematic 
Combination  of  these  two  effect  aggravates  (diagram  showing  the  formation  of 

the  discrepancy  between  the  prediction  and  ^  serrated  grain  boundary, 
the  experimental  observations. 


Lamellae  plate 
growing 
directioi 


The  analysis  above  indicates  that  the  step  height  is  another  very  important 
parameter  affecting  the  creep  resistance  for  fully  lamellar  spacing. 
Unfortunately,  this  parameter  is  barely  studied  in  the  literature  for  fully 
lamellar  y-TiAl  alloys,  and  not  available  for  either  of  the  four  materials  in 
Table  3-2.  Accordingly,  it  is  impossible  to  evaluate  its  effect  here 
quantitatively.  In  general,  the  step  height  may  depend  on  the  specific 
processing  route  used  in  the  preparation  of  the  materials  through  complex 
relationship.  It  is  evident  that  detailed  studies  on  this  aspect  are  necessary  to 
more  clearly  elucidate  the  experimental  data  in  terms  of  the  step  height  and 
the  interlamellar  spacing.  However,  it  is  worth  noting  that  the  step  height  in 
spray  formed  materials  may  have  more  chance  to  be  larger  than  that  in 
conventionally  processed  materials,  as  discussed  below. 


Qualitative  analysis  of  the  processing  effect  on  step  height  necessitates 
knowledge  on  the  fundamental  mechanisms  for  the  formation  of 
interlocking  at  the  neighboring  grain  boundaries.  No  existing  theory  is 
available  in  the  literature  regarding  the  formation  of  interlocking  for  fully 
lamellar  y-TiAl.  However,  studies  on  the  solidification  path  and  phase 
transformation  of  y-TiAl  indicated  that  a  lamellar  structure  is  formed  via 
either  of  the  following  reactions,  depending  on  the  composition  and  cooling 
rate  [12, 13, 18]: 

a  L(a  /  y)  ^  L(a2  /  y)  (low  composition,  high  cooling  rate) 


or 

a  — >  a2  L(a2  /  y)  (high  composition,  low  cooling  rate) 

where  L  stands  for  lamellar.  Similar  to  normal  single  phase  grain 
boundaries,  the  grain  boundaries  of  prior  a  (Ti  solid  solution,  bcc)  grains  or 
prior  a2  (TisAl,  hep)  grains  should  be  smooth.  Accordingly,  the  lamellae 
interlocking  could  only  evolve  during  precipitation  of  a2/y  lamellae  plates. 
It  is  reasonable  to  believe  that  the  evolution  of  the  lamellae  interlocking 
involves  encroach  of  the  earlier  precipitating  a2/y  lamellae  plates  (in  grain  1) 
into  later  precipitating  lamellae  plates  in  neighboring  grains  (grain  2),  as 
schematically  shown  in  Figure  3-12.  The  extent  of  encroachment,  or  the  step 
height,  is  therefore  determined  by  the  time  difference  between  the  initiations 
of  precipitation  of  these  two  lamellae  plates.  The  time  difference  may 
originate  from  the  variance  of  local  temperature  in  distinct  grains.  Compared 
with  the  conventional  processes,  variation  of  local  temperature  in  distinct 
grains  in  spray  forming  is  more  likely  to  occur,  since  in  spray  forming  the 
temperature  of  different  impinging  droplets  are  extremely  different  from  each 
other.  In  this  case,  serrated  steps  may  have  a  better  chance  to  evolve,  giving 
rise  to  larger  step  height. 


3.5.  Conclusions 

The  creep  behavior  of  spray  formed  y-TiAl  with  a  fine  equiaxed  fully 
lamellar  microstructure  was  studied  in  a  temperature-stress  regime  of  780  °C 
to  850  °C  and  180  MPa  to  320  MPa.  An  apparent  stress  exponent  of  4.3  and  an 
activation  energy  of  342  kj/mol  were  observed  in  the  high-temperature  high- 
stress  regime.  Compared  with  the  FL  y-TiAl  which  was  obtained  through 
conventional  casting  +  heat  treatment  processes,  the  spray  formed  y-TiAl 


exhibited  higher  creep  resistance.  The  higher  creep  resistance  observed  in  the 
present  study  was  discussed  in  light  of  the  interstitial  level,  the  chemical 
composition,  the  grain  size,  and  the  interlocking  of  lamellae  at  the  grain 
boundary  which  in  turn  may  be  a  function  of  interlamellar  spacing  and  the 
step  height  of  the  serrated  grain  boundaries.  It  was  proposed  that  the  small 
interlamellar  spacing  and  possibly  larger  step  height  may  contribute  to  the 
higher  creep  resistance  observed  in  the  present  study. 


CHAPTER  4. 


INVESTIGATION  OF  THERMAL  RESIDUAL  STRESSES  IN 
SPRAY  FORMED  y-TIAl  USING  FINITE  ELEMENT  METHOD 


4.1.  Objective 

The  purpose  of  this  study  is  to  investigate  the  thermal  residual  stresses  that 
develop  during  the  spray  deposition  processing  of  y-TiAl  using  thermo- 
elastoplastic  finite  element  analysis,  and  analyze  the  effect  of  different  cooling 
conditions  on  the  distribution  of  thermal  residual  stresses. 

4.2.  Finite  Element  Modeling 

With  the  increasing  availability  of  supercomputers,  the  finite  element 
method  is  becoming  an  increasingly  popular  approach  for  gathering  of 
information  which  can  not  be  extracted  from  experiments  alone,  in  a 
systematic  manner.  In  particular,  the  finite  element  method  has  generated 
considerable  interest  for  modeling  the  overall  thermo-elastoplastic  behavior 
of  materials  and  has  been  successfully  used  to  estimate  residual  stresses  in 
various  materials  [63-69].  The  theoretical  basis  for  the  thermo-elastoplastic 
finite  element  method  is  available  in  refs.  [70,  71]  in  detail. 

Since  the  stress  and  temperature  distributions  in  the  deposited  preform 
during  spray  atomization  and  deposition  process  are  interdependent,  the 
posed  thermo-mechanical  problem  is  fully  coupled.  Therefore,  the  thermal 
and  mechanical  solutions  must  be  performed  simultaneously.  In  the  present 
study,  ABAQUS,  a  commercially  available  finite  element  code  developed  by 
Hibbitt,  Karlsson,  and  Sorenson,  Inc.,  Pawtucket,  Rhode  Island,  was  employed 
due  to  its  capability  to  handle  a  coupled  temperature-displacement  problem. 
Elements  provided  in  ABAQUS  for  this  purpose  have  nodal  variables 
including  both  displacement  and  temperature.  These  elements  generally 
involve  a  low-order  interpolation  for  temperature  than  for  displacement 
(usually  parabolic  variation  of  displacement  and  linear  variation  of 
temperature)  in  order  to  obtain  a  compatible  variation  of  thermal  and 
mechanical  strains.  The  minimum  potential  energy  approach  to  the  finite 


element  formulation  of  displacement  problems  gives  the  following  general 
matrix  equation: 

F,  =  K^d  (4-1) 

where  F,  is  the  global  load  vector,  is  the  global  stiffness  matrix,  and  d  is 
the  global  displacement  vector.  The  application  of  Galerkin's  method  to  the 
thermal  problems  produces  the  following  matrix  equation: 

KJ  +  Ct  +  F,  =0  (4-2) 

where  T  is  the  temperature  vector,  is  the  conductivity  matrix,  C  is  the 
heat  capacity  matrix,  F^  is  the  heat  flow  vector.  The  heat  capacity, 
conductivity,  and  stiffness  matrices  in  equations  (4-1)  and  (4-2)  are  generated 
using  a  full  Gaussian  integration  scheme,  and  the  coupled  system  is  solved 
using  Newton's  method  [72].  The  detailed  description  of  solution  procedures 
for  a  coupled  temperature-displacement  analysis  is  available  in  reference  [73]. 

4.2.1.  Materials 

In  the  present  study,  gamma  titanium  aluminide  alloy  (y-TiAl)  based  on  Ti- 
47  at.%  A1  was  selected  as  a  model  material.  The  phase  diagram  of  titanium- 
aluminum  is  shown  in  Figure  4-1  [74].  Gamma  TiAl  has  special  interest  in 
aerospace  and  automotive 
industries  due  to  its  low  density, 
high  melting  temperature, 
elevated  temperature  strength, 
high  resistance  to  oxidation  and 
hydrogen  absorption,  and  excellent 
creep  properties  [12,  13,  16,  75,  76]. 

However,  y-TiAl  suffers  from  low 
ductility  and  toughness  at  ambient 
temperatures  which,  along  with 
poor  formability,  appears  to  be  the 
single  most  serious  obstacle  to  its  fabrication  and  full  utilization  [77].  In  order 
to  improve  the  ambient  temperature  ductility  of  y-TiAl,  powder  metallurgy 
or  rapid  solidification  processes  followed  by  powder  consolidation,  an 
approach  which  promotes  a  grain  refinement,  is  typically  employed. 
Considering  the  fact  that  the  multiple  steps  of  powder  production,  sieving, 
de-gassing,  and  consolidation  may  be  consolidated  into  a  single  processing 
step,  spray  atomization  and  deposition  process  stands  out  as  a  potentially 
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Figure  4-1.  Phase  diagram  of  titanium- 
aluminum  [74]. 


economical  alternative  to  produce  ductility-improved  y-TiAl.  The  physical 
and  mechanical  properties  of  y-TiAl  at  various  temperatures  used  in  the 
present  study  are  listed  in  Table  4-1  [75,  78-81]. 


4,2.2.  Geometry 

In  the  present  investigation,  due  to  the  difficulty  in  modeling  the  transient 
problem  posed  by  a  growing  preform  during  spray  atomization  and 
deposition,  the  fully  grown  preform  at  the  end  of  deposition  was  modeled. 
Accordingly,  the  present  model  formulation  {i.e.,  initial  conditions)  begins 
immediately  following  the  completion  of  the  deposition  stage.  The 
geometrical  shape  of  the  final  preform  considered  in  this  study  was  assumed 
to  exhibit  a  contour  akin  to  the  Gaussian  distribution  of  droplets  impacting 
on  the  substrate  as  shown  in  Figure  4-2.  The  height  and  the  bottom  diameter 
of  the  deposited  preform  were  assumed  to  be  27  and  100  mm,  respectively, 
consistent  with  available  experimental  results  [82,  83]. 


Table  4-1.  The  physical  and  mechanical  properties 
of  y-TiAl  at  various  temperatures  [75,  78-81] 


Temperature 

(K) 

Elastic 

Modulus 

(GPa) 

Yield 

Stress 

(MPa) 

Specific 

Heat 

(J/kg-K) 

Thermal 

Expansion 

(prn/m) 

Thermal 

Conductivity 

(W/m-K) 

300 

170 

410 

658 

12.87 

89.86 

473 

168 

400 

738 

13.77 

89.86 

673 

160 

390 

780 

14.40 

89.86 

873 

154 

350 

808 

15.00 

89.86 

1073 

147 

300 

807 

15.60 

37.57 

1273 

138 

150 

805 

16.20 

37.57 

1473 

129 

75 

800 

16.80 

37.57 

1753 

115 

38 

800 

17.50 

37.57 

Density  (kg/m? ) :  3852.5 

As  shown  in  Figure  4-2,  the  deposited  preform  exhibits  an  axisymmetric 
shape.  Therefore,  the  three-dimensional  preform  could  be  modeled  using 
two-dimensional  axisymmetric  elements.  Initially,  a  coarse  mesh  using 
relatively  small  number  of  elements  was  employed  for  modeling,  and  the 
number  of  elements  was  gradually  increased  to  obtain  convergent  solutions. 
A  total  of  486  four-noded  temperature-displacement  elements  was  used  in 


this  study  [72].  Figure  4-3  shows  the  two-dimensional  axisymmetric  element 
mesh  used  herein  for  modeling  the  preform. 

Top  Layer 


Figure  4-2.  The  geometrical  shape 
of  deposited  preform  from  spray 
atomization  and  deposition. 
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z 

Figure  4-3.  Two-dimensional 
axisymmetric  element  mesh  for 
modeling  of  preform. 


4.2.3.  Boundary  and  initial  conditions 

Throughout  the  modeling,  symmetric  boundary  conditions  were  imposed 
such  that  every  symmetric  line  is  constrained  in  its  normal  direction  to 
represent  the  effect  of  the  rest  of  the  material.  For  the  initial  temperature 
condition  of  the  modeling,  the  temperatures  at  the  top  and  bottom  surfaces  of 
the  preform  were  assumed  to  be  1753  and  1250  K,  respectively,  since  the  top 
surface  of  the  preform  was  covered  with  partially  solidified  droplets  with  a 
solidus  temperature  of  1753  K  at  the  end  of  the  deposition  process,  and  the 
bottom  surface  of  the  preform  was  cooled  down  during  the  process  and 
assumed  to  have  a  temperature  of  1250  K,  when  the  deposition  process  was 
completed.  The  temperature  inside  the  preform  was  then  linearly 
interpolated. 

During  the  spray  atomization  and  deposition  process,  the  deposited  preform 
loses  thermal  energy  through  a  combination  of  three  heat  transfer  processes: 
a)  conduction  to  the  low  temperature  substrate,  b)  convection  to  the 
atomization  gas,  c)  radiation  to  the  environment  [83].  On  the  basis  of  thermal 
flux  measurements,  Mathur  et  al.  [84]  concluded  that  the  heat  transfer  process 
during  the  initial  stage  of  deposition  was  very  high  and  that  the 
corresponding  thermal  conductive  heat  transfer  coefficient  between  deposited 
preform  and  substrate  was  approximately  10^  W/m^K.  Following  a  brief 
transient  period,  the  magnitude  of  the  heat  transfer  coefficient  decreases  to 
approximately  550  W/m^K  under  steady  state  condition  for  the  copper 
substrate  cooled  by  water.  The  convective  heat  transfer  coefficient  from  the 


top  surface  of  the  preform  to  the  atmosphere  at  the  end  of  deposition  is 
assumed  to  be  at  least  one  order  of  magnitude  smaller  than  the  conductive 
heat  transfer  coefficient  from  the  bottom  surface  of  the  preform  to  the 
substrate  [27].  Therefore,  in  this  study,  heat  transfer  coefficients  of  21  and  550 
W  /  m  2K  were  used  for  the  top  and  bottom  surfaces  of  the  preform, 
respectively,  when  the  substrate  material  was  copper,  which  represents  the 
highly  conductive  substrate  case.  In  addition,  in  order  to  investigate  the  effect 
of  different  cooling  conditions  on  the  thermal  residual  stresses,  the  same 
analysis  was  performed  using  two  different  conductive  heat  transfer 
coefficients,  which  represent  for  an  insulated  substrate  case  and  an 
intermediate  case,  respectively,  with  the  same  convective  heat  transfer 
coefficient  of  21  W/m^K.  The  heat  transfer  coefficients  used  in  this  study  are 
listed  in  Table  4-2.  The  radiative  heat  transfer  was  not  considered  in  this 
study,  since  its  effect  becomes  significant  only  when  the  temperature  of  the 
deposited  preform  is  extremely  high  [83].  Throughout  the  modeling  of  three 
different  cases,  the  same  initial  temperature  condition  of  the  preform  was 
used  in  order  to  focus  on  the  effect  of  the  substrate  with  different  thermal 
conductivities.  The  ambient  temperature  was  assumed  to  be  300  K. 


Table  4-2.  Heat  transfer  coefficients  used  for  the  modeling 


Substrate  Material 

Conductive  Heat  Transfer 
Coefficient  (W/m^K) 

Convective  Heat  Transfer 
Coefficient  (W/m^K) 

Highly  Conductive  (Cu) 

550 

21 

Intermediate 

275 

21 

Insulated 

0 

21 

4.3.  Results  and  Discussion 

Figures  4-4,  4-5,  and  4-6  show  the  von  Mises'  equivalent  stresses  along  the 
top,  middle,  and  bottom  layers  in  the  preform  as  indicated  in  Figure  4-3.  For 
the  conductive  substrate  cases,  the  equivalent  stresses  have  higher  values  for 
the  region  closer  to  the  outer  edge  of  the  preform.  In  other  words,  the  stresses 
increase  as  the  deposit  thickness  decreases.  This  is  due  to  the  fact  that  the 
region  with  lower  deposit  thickness  experiences  higher  thermal  gradient 
during  cool  down.  The  region  closer  to  the  edge  of  the  preform  has  stresses 
that  are  greater  than  410  MPa,  which  is  the  yield  stress  of  y-TiAl  at  300  K. 


Therefore,  plastic  deformation  or  possibly  cracks  may  develop  in  that  region. 
For  the  insulated  substrate  case,  however,  the  equivalent  stress  distribution  is 
significantly  reduced  and  exhibits  a  relatively  homogeneous  pattern  for  the 
three  different  layers  as  compared  with  those  of  two  other  cases.  Overall,  it 
was  found  that  the  equivalent  stress  distribution  is  a  strong  function  of 
cooling  rate  in  the  preform  and  the  deposit  thickness.  It  was  also  found  that 
the  equivalent  stresses  decrease  as  the  conductive  heat  transfer  coefficient  (h) 
decreases. 
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along  the  top  layer  in  the  preform  for 
three  different  cooling  conditions. 
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Figure  4-5.  The  von  Mises'  stresses 
along  the  middle  layer  in  the  preform  for 
three  different  cooling  conditions. 


Figure  4-7  shows  the  equivalent  stresses  at  the  center  of  the  preform  as  a 
function  of  distance  along  the  z  direction  from  the  top  to  bottom  of  the 
preform.  In  the  center  region,  it  was  found  that  the  stress  reaches  the 
maximum  value,  and  the  stress  distribution  decreases  as  the  h  decreases. 
Depending  on  the  h  value,  the  maximum  stress  varies  from  110  to  280  MPa. 
It  was  also  noticeable  that  as  the  h  decreases,  the  stresses  in  the  center  region 
gradually  tend  to  be  flattened. 

The  radial,  hoop,  and  axial  stresses  along  the  center  line  of  the  preform  are 
shown  in  Figures  4-8,  4-9,  and  4-10,  respectively.  The  radial  stresses  are 
mostly  tensile  for  the  upper  region  of  the  preform  and  compressive  for  the 
lower  region  of  the  preform  for  all  three  cases  considered  in  this  study.  A 
similar  trend  can  be  observed  for  the  hoop  stresses.  The  axial  stresses, 
however,  are  compressive  for  the  entire  preform.  The  magnitude  of  these 


stresses  increases  as  the  distance  from  the  top  of  the  preform  increases  {i.e., 
closer  to  the  bottom  of  the  preform).  Again,  the  radial,  hoop,  and  axial 
stresses  decrease  as  the  thermal  conductivity  of  the  substrate  decreases. 
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Figure  4-6.  The  von  Mises'  stresses 
along  the  bottom  layer  in  the  preform  for 
three  different  cooling  conditions. 
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Figure  4-7.  The  von  Mises'  stresses  at 
the  center  of  the  preform  as  a  function  of 
distance  along  the  Z  direction  from  the 
top  to  bottom  of  the  preform  for  three 
different  cooling  conditions. 
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Figure  4-8.  The  radial  stresses  at  the 
center  of  the  preform  as  a  function  of 
distance  along  the  Z  direction  from  the 
top  to  bottom  of  the  preform  for  three 
different  cooling  conditions. 
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Figure  4-9.  The  hoop  stresses  at  the 
center  of  the  preform  as  a  function  of 
distance  along  the  Z  direction  from  the 
top  to  bottom  of  the  preform  for  three 
different  cooling  conditions. 


Figures  4-11,  4-12,  and  4-13  show  the  radial,  hoop,  and  axial  stresses  along  the 
middle  layer  of  the  preform.  The  radial  stresses  in  this  region  are  mostly 
tensile  for  the  entire  preform,  and  these  values  are  greatly  reduced  when  the 
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Figure  4-10.  The  axial  stresses  at  the 
center  of  the  preform  as  a  function  of 
distance  along  the  Z  direction  from  the 
top  to  bottom  of  the  preform  for  three 
different  cooling  conditions. 
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Figure  4-11.  The  radial  stresses  along 
the  middle  layer  of  the  preform  for  three 
different  cooling  conditions. 
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Figure  4-12.  The  hoop  stresses  along 
the  middle  layer  of  the  preform  for  three 
different  cooling  conditions. 
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Figure  4-13.  The  axial  stresses  along 
the  middle  layer  of  the  preform  for  three 
different  cooling  conditions. 


insulated  substrate  case  is  considered.  On  the  other  hand,  the  hoop  stresses 
are  tensile  for  the  inner  area  of  the  preform  and  compressive  for  the  outer 
area  of  the  preform  for  the  conductive  substrate  cases.  When  the  substrate  is 
insulated,  the  hoop  stresses  become  tensile  for  the  entire  preform.  Contrary 
to  the  hoop  stress  distribution,  the  axial  stress  distribution  exhibits 
compressive  values  for  the  inner  area  of  the  preform  and  tensile  values  for 
the  outer  area  of  the  preform. 


4.4.  Summary 

Based  on  the  present  study,  the  following  conclusions  were  made: 

1.  Thermal  residual  stresses  were  developed  in  the  spray  formed  y-TiAl 
during  cooling  from  the  deposition  temperature  to  ambient 
temperature. 

2.  The  axial,  hoop,  and  radial  stresses  were  found  to  be  mostly 
compressive  for  the  lower  center  region  of  the  preform.  The 
magnitude  of  these  compressive  stresses  increases  for  the  region  closer 
to  the  bottom  of  the  preform.  For  the  other  region,  the  axial,  hoop,  and 
radial  stresses  become  tensile. 

3.  The  residual  stress  distribution  is  a  strong  function  of  cooling  rate  in 
the  deposited  preform,  and  the  residual  stresses  decrease  as  the  cooling 
rate  decreases.  When  the  substrate  is  insulated,  the  preform  has  the 
lowest  residual  stress  distribution  for  the  cases  studied  herein. 

4.  The  residual  stress  distribution  is  a  strong  function  of  deposit 
thickness.  For  the  region  closer  to  the  outer  edge  of  the  preform,  which 
experiences  the  most  severe  thermal  gradient,  has  the  highest  residual 
stresses.  The  von  Mises'  equivalent  stresses  in  this  region  were  found 
to  be  greater  than  the  yield  stress  of  y-TiAl.  Therefore,  it  is  anticipated 
that  plastic  deformation  or  cracks  may  be  present  in  this  region  of  the 
preform. 


CHAPTER  5. 


TWO  DIMENSIONAL  MODELING  OF  MOMENTUM  AND  THERMAL 
BEHAVIOR  DURING  SPRAY  ATOMIZATION  OF  y-TiAl 


5.1.  Objective 

The  present  study  was  undertaken  with  the  objective  of  formulating  a  two- 
dimensional  model  to  investigate  the  thermal  and  momentum  behavior  of 
y-TiAl  during  spray  atomization.  A  specific  alloy  composition  of  Ti-48  at.%Al 
was  selected  for  study  due  to  recent  development  of  commercial  alloys  based 
on  this  stoichiometry:  Ti-48Al-2.5Nb-0.3Ta  [2],  Ti-48Al-2Nb-2Zr  [2],  Ti-48A1-2V 
[14]  and  Ti-48Al-2Cr-2Nb  [85]  (All  are  in  atomic  percentage). 


5.2.  Model  Formulation 


The  primary  assumptions  underlying  the  problem  formulation  that  will  be 
discussed  in  subsequent  sections  are  listed  as  follows: 

(1) .  The  alloy  melt  is  atomized  instantaneously  in  a  small  atomization  zone, 

whose  magnitude  is  approximately  3~4  times  the  inside  diameter  (~0.001 
m)  of  a  metal  delivery  tube. 

(2) .  Following  atomization,  the  concentration  of  the  droplets  in  the  gas 

environment  is  sufficiently  dilute,  so  that  the  droplet-droplet 
interactions  may  be  safely  neglected. 

(3) .  The  flow  of  the  atomization 

gas  is  axisymmetric.  As  a 
result,  the  droplets  do  not 
experience  rotation  around 
the  spray  cone  axis.  The 
motion  of  the  droplets  is  then 
assumed  to  be  two 
dimensional.  Accordingly,  a 
radial  coordinate  system,  (r, 

Z),  is  selected  for  the  present 

problem,  where  r  represents  the  radial  distance  and  Z  represents  the  axial 
distance.  The  coordinate  system  is  shown  schematically  in  Figure  5-1. 


Figure  5-1.  Schematic  representation  of  the 
spray  atomization  and  deposition  process. 
The  radial  coordinate  system,  (r,  Z),  is 
shown  in  the  insert  on  the  right  hand  side. 


(4).  The  (r,  Z)  coordinates  at  the  center  of  the  atomization  zone  are  selected 
to  correspond  to  (Om,  Om),  where  the  radial  droplet  velocity  is  zero,  and 
the  axial  droplet  velocity  has  a  relatively  small  value  (for  example,  1 
m/s)  as  a  result  of  the  initial  free  fall  of  the  alloy  melt.  Consistent  with 
assumption  (1),  the  initial  positions,  (ro,  Zq),  of  atomized  droplets  will  be 
in  the  range  from  (Om,  Om)  to  (Om,  0.004m).  In  addition,  it  is  also 
assumed  in  the  present  study  that  droplets  are  randomly  distributed  in 
the  atomization  zone  once  they  are  formed. 

5.2.1.  Momentum  profile  of  atomization  gas 

Based  on  the  experimental  results  by  Bewlay  et  al.  [86]  and  Grant  et  al.  [87],  in 
which  the  geometry  of  the  close-coupled  atomization  device  used  is  similar  to 
that  used  in  the  present  investigation,  the  velocity  profile  of  the  atomization 


gas  may  be  described  as  follows  [88, 89]: 

u  2g  (0,Z)  =  u  gQ  exp  (  -  Z/  X)  (5-1) 

with  ^  =  3.04x10  (5-2) 

gO 

and  t)2g(r,Z)  =  'i)2g(0,Z)[l  —  (q  qq4  ^  0.268Z^  ^  "^^oo  (5-3) 

rg('^’^)  ~  ^  zg  (*'’Z)(  Q  Q  +  Z  ^  (5-4) 


where  Uzg  and  Urg  are  the  axial  and  radial  components  of  gas  velocity, 
respectively;  Ugo  is  the  initial  gas  velocity;  is  the  axial  gas  velocity  at  Z=  oo; 
and  X  is  the  exponential  decay  coefficient. 

Eqns.  (5-3)  and  (5-4)  are  applicable  only  for  r/(0.004+0.268Z)<l,  i.e.,  inside  the 
spray  cone.  When  r/(0.004+0.268Z)>l,  it  is  reasonable  to  assume  that  \)zg(r,  Z) 
and  'Urg(r,  Z)  are  both  close  to  zero. 

5.2.2.  Velocity  profile  of  droplets 

During  atomization,  droplets  are  subjected  to  an  accelerating  or  decelerating 
drag  force  as  a  result  of  the  differences  in  velocity  between  droplets  and 

atomization  gas  [90]: 
dt) 

=  V(P,-  p,)g  +  (A/8)C„„P.(1),.  (5-5) 


(5-6) 


where  pd  and  Pg  are  the  densities  of  droplet  and  atomization  gas  (Ar), 
respectively;  V  and  A  are  the  volume  and  surface  area  of  droplet,  respectively; 
Uzd  and  ^rd  are  the  axial  and  radial  components  of  droplet  velocity, 
respectively;  Udg  (=[('Ozg-nzd)^+(nrg-nrd)^]^/^)  is  the  relative  velocity  between 
gas  and  droplet;  Cdrag  is  the  drag  coefficient;  and  g  is  the  gravity  constant. 

Droplet-droplet  interactions.  Basset  history  term  and  added  mass  effect  [90]  are 
neglected  in  Eqns.  (5-5)  and  (5-6)  since  these  terms  are  generally  much  smaller 
than  the  inertia  and  the  drag  forces  in  processes  such  as  spray  atomization. 

A  well  known  relationship  between  the  drag  coefficient,  Cdrag/  arid  the 
Reynolds  number,  Nre  (=DpgUdg/Pg,  where  D  is  the  droplet  diameter  and  |ig 
the  gas  viscosity),  is  =  0.28  +  (6/N^^)  +  (21/Nj^^)  [91],  which  is  applicable 

for  0.1<NRe  <4000.  In  the  present  study,  Nrc  falls  into  a  wider  range  (i.e. 
0.02 — 8000)  than  the  one  mentioned  above.  Therefore,  an  improved 
expression  for  Cdrag/  as  a  function  of  Nrb,  may  be  used  as  described  in 
reference  [89]: 


r24/N„  . 

. 0<N^  <1 

Re 

Re 

. 

. 1<N„  <400 

Re 

(5-7) 

0.5 . 

. 400  <N^  <3x  10^ 

Ke 

Using  Eqns.  (5-5)-(5-7),  together  with  the  gas  velocity  profiles  calculated  from 
Eqns.  (5-l)-(5-4)  and  the  parameters  listed  in  Table  5-1,  the  two-dimensional 
trajectory  and  velocity  profile  of  a  given  droplet  may  be  calculated 
numerically  once  ro  and  'Uzd(r/0)  are  given. 

5.2.3.  Thermal  and  solidification  history  of  droplets 

In  terms  of  Tj,  Tn,  Tl  and  Tp,  which  represent  initial,  nucleation,  liquidus  and 
peritectic  temperatures  of  y-TiAl,  respectively,  the  thermal  and  solidification 
profiles  of  y-TiAl  may  be  classified  into  five  well  defined  stages  (refer  to 
Figure  5-2  and  Figure  5-7):  (1)  cooling  in  the  liquid  state  when  Tn<T<Ti,  (2) 
recalescence  when  T=Tn,  (3)  segregated  solidification  when  Tp<T<Tl/  (4) 
peritectic  solidification  when  T=Tp,  and  (5)  cooling  in  the  solid  state  following 
the  completion  of  solidification.  Before  any  detailed  calculation  of  the  droplet 


thermal  and  solidification  histories  can  be  accomplished,  the  values  of  Tj,  Tn, 
Tl  and  Tp  need  to  be  determined  first.  Tj  may  be  determined  from  the  degree 
of  melt  superheat;  Tp  and  Tp  may  also  be  easily  obtained  according  to  Ti-Al 
phase  diagram  (Figure  5-2).  The  determination  of  Tn,  however,  requires 
detailed  knowledge  of  the  relevant  solidification  kinetics.  Relevant 
formulation  is  described  in  more  detail  below. 


Table  5-1.  Physical  parameters  used  in  Eqns.  (5-1)  through  (5-20) 


parameter 

value 

source 

^pg 

518.8  Jkg-iK-i 

[92,93] 

Cpi 

829  Jkg-iK-i 

RMA’' 

Cps 

727  Jkg-iK-i 

RMA 

Hf 

28.27  kjmol-i 

[Eqn.(5-21)] 

ke 

0.74 

[Kg.  2] 

0.01636  Wm-iK-i 

[92,93] 

R 

0.01  ms-iK-i 

[94] 

Ti 

1885  K 

[12] 

Tl 

1785  K 

[12] 

Tm 

2016  K 

[12] 

Tp 

1764  K 

[12] 

T 

Oi 

273  K 

[92] 

T 

^  Of 

473  K 

[84] 

a 

0.1  m 

[88] 

DgO 

150  ms'i 

[86,88] 

0.004  ms"' 

[86] 

8 

0.5 

[89] 

1.7  X  lO'S  kgm‘is‘1 

[92] 

ftm 

3.6  X 10-3  kgm'is'i 

RMA 

Pd 

3636  kgm-3 

RMA 

Pg 

1.654  kgm'3 

[92,93] 

O'!!! 

203.4  X 10-6  kJm-2 

[Eqn.  (5-22)] 

a 

2.08  X 10-6  m^mol'i 

RMA 

*Rule  of  mixtures  approach. 


Based  on  the  work  of  Hirth  [95]  and  Turnbull  [96],  the  maximum 
homogeneous  undercooling  AThom  (=Tl"Tn)  for  a  given  droplet  may  be 
calculated  as  follows: 


(5-8) 


^  ^  -L< 

■  3k(T  ,  -  AT,^  [IO^VaT^^^/T] 

where  Hf,  Om  and  ^  are  the  latent  heat  of  fusion,  the  solid-liquid  interface 
energy  and  the  atomic  volume  of  y-TiAl,  respectively;  k  is  the  Boltzmann's 
constant;  V  is  the  droplet  volume  in  m^;  and  t  is  the  cooling  rate  of  droplet 
in  the  first  stage,  namely  cooling  in  the  liquid  state. 


As  may  be  seen  later  (refer  to  Section  5.4.4),  T  depends  strongly  on  the  droplet 
size  (D),  initial  position  (ro),  and  axial  distance  from  the  atomization  zone  (Z). 
As  an  approximation,  a  value  averaged  on  different  ro  and  Z  was  used  for 
each  droplet  considered  in  present  study.  This  is  justified  by  results  of 
numerical  solution  of  Eqn.  (5-8),  which  indicate  that  the  value  of  AThom  is 


relatively  insensitive  to  the  changes  in 


Figure  5-2.  Central  part  of  the  equilibrium 
Ti-Al  phase  diagram  [12,  18].  Dotted  lines 
are  drawn  for  calculation  of  ke  and  Tm  in 
equations  (5-18)  and  (5-19).  Dashed  line 
represents  the  alloy  composition  in  the 
present  study. 


T. 


Figure  5-3.  Relationship  between  actual 
undercooling,  homogeneous  undercooling 
and  droplet  size. 


It  is  well  established,  nevertheless,  that  heterogeneous  rather  than 
homogeneous  nucleation  mechanisms  limit  the  degree  of  undercooling 
experienced  by  a  large  proportion  of  the  droplet  population  during 
atomization  [84,  97].  Only  in  very  small  droplets  (e.g.,  D<20  |xm)  does 
homogeneous  nucleation  play  an  important  role  during  solidification. 
Accordingly,  the  actual  degree  of  undercooling  experienced  by  a  droplet, 
ATactuai/  will  generally  be  significantly  smaller  than  the  undercooling 
necessary  for  homogeneous  nucleation,  AThom-  Based  on  experimental  results 


of  Mathur  et  al.  [84],  an  exponential  relationship  between  ATactuai  and  AThom 

may  be  formulated  as  (Figure  5-3): 

ATactuai  =  AThomexp(-2.2  X 1012V)  (5-9) 

It  should  be  noted  that  the  degree  of  undercooling  that  is  present  during 
heterogeneous  nucleation  is  controlled  by  the  availability  of  heterogeneous 
nucleation  sites  in  the  volume  or  surface  of  the  droplet  and  the  wetting  angle 
between  nucleant  and  melt.  However,  the  availability  of  heterogeneous 
nucleants  is  more  important  than  the  wetting  angle  between  nucleant  and 
melt.  The  former  strongly  depends  on  the  droplet  size  [98,  99].  Accordingly, 
ATactuai/AThom  for  a  specific  droplet  during  atomization  is  more  significantly 
affected  by  droplet  size  than  by  alloy  composition  (AThom  has  already  taken 
into  account  some  of  the  effects  derived  from  the  alloy  composition). 
Actually,  the  experimental  results  obtained  by  Mathur  et  al.  [84]  reveal  that 
ATactuai/AThom  is  relatively  insensitive  to  chemical  composition  of  the 
materials  sprayed  (see  Figure  5-3).  Therefore,  the  applicability  of  Eqn.  (5-9)  is 
extended  to  the  y-TiAl  studied  herein,  despite  the  fact  that  this  correlation 
was  originally  formulated  on  the  basis  of  experimental  results  obtained  for 
Ni-20  wt.%  Cr  and  Fe-20  wt.%  Mn. 

Eqns.  (5-8)  and  (5-9)  may  be  used  to  calculate  ATactuai  arid  then  Tn  (=Tl- 
ATactuai)  for  arty  droplet  with  a  given  size  (D  or  V).  The  values  of  ATactuai 
calculated  for  the  range  of  droplet  sizes  studied  herein  are  listed  in  Table  5-2. 
Generally,  experimental  measurement  of  the  magnitude  of  undercooling 
experienced  by  droplets  during  spray  atomization  is  very  difficult.  To  that 
effect,  inspection  of  available  literature  indicates  that  there  is  no  experimental 
data  reported  so  far  on  the  undercooling  experienced  by  atomized  y-TiAl 
droplets  during  spray  atomization.  It  is  worth  noting,  however,  that  in  a 
related  study  conducted  by  Valencia  et  al.  [100],  an  undercooling  as  large  as  286 
K  was  experimentally  achieved  using  a  levitation  technique  for  a  Ti-50A1 
droplet  with  a  diameter  of  approximately  5000  pm. 


Table  5-2.  Calculated  actual  undercooling 
of  y-TiAl  with  various  droplet  sizes 


size(pm) 

20 

50 

80 

175 

325 

500 

undercooling(K) 

469 

393 

247 

1 

5.2. 3.1.  Cooling  in  the  liquid  state  (Tm<T<Ti) 

In  a  completely  molten  droplet,  thermal  energy  is  dissipated  into  the 
environment  via  a  combination  of  convection  and  radiation  until  Tn  is 
reached: 

-  Pd''Cp,^=  h(T  -  T„)A  +  gc<t‘'  -  ry  (5-10) 

where  Cpi  is  the  specific  heat  of  liquid  y-TiAl;  £  is  the  emissivity  of  the 
surface;  a  is  the  Stephen-Boltzmann  constant;  h  is  the  convection  heat 
transfer  coefficient,  which  may  be  expressed  as  [101]: 

h  =  (kg/D)[2  +  0.6N|;^^Np"]  (5-11) 

where  Npr  (=  |lgCpg/kg)  is  the  Prandtl  number;  kg  and  Cpg  are  the  thermal  conductivity  and 
the  specific  heat  of  gas,  respectively;  and  To  is  the  temperature  of  the  gas,  which  is 
assumed  to  increase  exponentially  as  a  result  of  heat  transfer  between  the  gas  and  the 
droplets  [88]: 

To=T„t-(To,-T„i)exp(-Z/a)  (5-12) 

where  Tq;  and  Tof  are  initial  and  final  gas  temperatures,  respectively;  and  a  is  a  constant. 

The  approximation  of  Newtonian  cooling  underlying  Eqn  (5-10)  is  acceptable 
when  the  magnitude  of  the  Biot  number,  Bi  (=  [h(D/2)]/kd,  where  k^  is  the 
thermal  conductivity  of  droplet),  is  less  than  0.1  [92, 102].  In  the  present  study, 
Bi<0.012. 

Eqns.  (5-10)-(5-12),  along  with  the  values  of  Tj  and  Tn  may  be  used  to  predict 
the  temperature,  T,  for  a  given  droplet  at  any  time,  t,  or  any  position,  (r,  Z).  In 
this  stage  the  fraction  of  solid,  f,  contained  in  a  droplet  is,  by  definition,  zero. 

5.2. 3.2.  Recalescence  (T=Tm) 

Recalescence  is  the  rapid  heating  of  an  undercooled  molten  sample  due  to  the 
nucleation  of  the  solid  phase.  Since  recalescence  is  driven  by  the  magnitude 
of  the  undercooling  present,  as  the  droplet  heats  up  the  rate  of  heat  release 
decreases  and  recalescence  terminates  when  the  rate  of  heat  release  becomes 
equal  to  the  rate  of  heat  extraction  from  the  outer  surface  of  the  droplet: 

Hj[p_,V(l-t,)l^  =  h(T,-T^)  (5-13) 

where  fr  and  Tr  are  the  fraction  solidified  in  the  droplet  and  the  temperature 
of  the  droplet  at  the  end  of  recalescence,  respectively. 


The  rate  of  growth  of  solid  phase  in  undercooled  melt  with  small  to  moderate 
undercoolings  may  be  given  by  [94]: 


df  _  RAT 
dt  D 


(5-14) 


in  which  R  is  the  solid/liquid  interface  mobility,  and  AT=TL-Tr  at  the  end  of 
recalescence. 


By  establishing  a  thermal  energy  balance  between  the  latent  heat  released  and 

the  thermal  energy  that  is  needed  to  accommodate  the  temperature  rise  in  the 

droplet,  fr  in  Eqn.  (5-13)  may  be  calculated  as: 

Hffr  =  (Tr  -  TN)[Cpi(l-fr)  +  Cpsfr]  (5-15) 

where  Cps  is  the  specific  heat  of  solid  y-TiAl. 

Using  Eqns.  (5-13)-(5-15),  the  droplet  temperature  and  the  fraction  solidified  at 
the  end  of  recalescence  may  be  readily  determined.  Detailed  calculations 
indicate  that,  when  (h/(pdCpsR))<2  K,  Eqns.  (5-13)-(5-15)  may  be  approximated 
as 

T;-{-^  +  T^+T^)T^  +  -^T^  +  TJ^  =  0  (5-16) 

^  ps 


Taking  into  account  the  fact  that  Tj.  <  Tj^,  then  Eqn.  (5-16)  is  equivalent  to 

(5-17) 


T  = 

^  r 


Tl- . (H/Cp,)4.T^>T^ 


T„  +  (Hj/C 


ps 


). 


•(H/Cp3)-hT^<Tl 


In  the  present  study,  the  values  of  h/(pdCpsR)  are  smaller  than  1.41  K  for 
droplet  sizes  corresponding  to  D>20  |im. 

In  addition,  available  results  reveal  that  the  time  interval  that  is  associated 
with  the  recalescence  phenomenon  is  less  than  seconds  [103].  Therefore, 
it  is  safe  to  assume  that  the  recalescence  is  essentially  complete  within  a  single 
calculation  step  in  present  study  (which  is  in  the  order  of  milliseconds). 

5.2. 3.3.  Segregated  solidification  (Tp<T^i_) 

The  segregation  behavior  of  the  solute  following  the  completion  of 
recalescence  may  be  described  by  the  Scheil  equation  [104],  the  Brody-Flemings 
theory  [105],  the  Bower-Brody-Flemings  theory  [106]  and  the  Clyne-Kurz 


theory  [107].  However,  it  was  found  that  the  selection  of  a  particular  theory 

has  only  a  limited  effect  on  the  overall  cooling  behavior  of  atomized  droplets 

[88].  Therefore,  because  of  its  simplicity,  the  Scheil  equation  is  employed  to 

describe  segregation  behavior  in  the  present  study: 

(ci-cs)df  =  (l-f)dci  (5-18a) 

or,  after  being  integrated, 

c,-30  l/(k-l) 

(5-18b) 

where  cq,  Cg  and  ci  are  the  initial ,  solid  and  liquid  composition  of  solidifying 
y-TiAl  droplet,  respectively;  and  kg  [=(cs-30)/(ci-30),  refer  to  Figure  5-2]  is  a 
constant  similar  to  the  conventional  equilibrium  coefficient. 


Writing  ci  and  cq  in  terms  of  T,  Tl  and  Tm  yields: 

T  -T  i/(k,-i) 


M 


and 


df 

dT 


1-f, 


-G 


M 


2-k 

T 

-) 


(ke-l)(TM-TL)^T^-T^ 

where  Tm  is  a  reference  temperature  as  defined  in  Figure  5-2. 


(5-19a) 

(5-19b) 


The  thermal  energy  balance  governing  the  thermal  history  of  the  droplet 
during  this  stage  is  identical  to  that  described  by  Eqn.  (5-10)  if  Cpi  is  substituted 

An 

with  C  ,(1  -  f)  +  in  which  the  third  term  takes  into  account  the 

pi^  '  ps  fdT 

latent  heat  released  during  solidification. 


5.23.4.  Peritectic  solidification  (T=Td) 

In  this  stage,  the  droplet  temperature  remains  at  the  peritectic  temperature, 
Tp,  until  the  droplet  is  completely  solidified.  It  then  follows  that  the  fraction 
of  solid  present  in  the  droplet,  as  a  function  of  time  or  position  may  be 
expressed  as: 

PdH  =  h(T  p  -  Tq)  A  ea(Tp  -  T  A  (5-20) 


5.2. 3.5.  Cooling  in  the  solid  state  (f-1) 

In  this  stage,  f=l  and  the  associated  thermal  energy  balance  is  established  by 
using  Eqn.  (5-10),  and  substituting  Cpg  for  Cpi. 


5.3.  Numerical  Solution 


Computer  programs  have  been  prepared  for  the  execution  of  the  calculations 
required  by  Eqns.  (5-1)  through  (5-20).  The  primary  point  of  the  present 
computational  exercise  is  to  concentrate  on  a  given  droplet,  formed  near  the 
atomization  point  (e.g.  Zo=0  m,  ro=0.001  m),  having  a  particular  diameter  (say 
80  pm),  and  then  follow  its  two-dimensional  trajectory  as  it  travels  down  the 
atomization  chamber. 

Numerical  calculation  will  be  carried  out  for  each  pair  of  droplet  size  and 
droplet  initial  position,  {D,  (Zq,  ro)}.  Consistent  with  the  results  of  Moll  et  al 
[108],  the  characteristic  sizes,  Di6,  D50  and  034,  which  represent  droplet 
diameters  for  which  16  wt.  %,  50  wt.  %  and  84  wt.  %  of  droplets  are  smaller 
than  these  values,  respectively,  were  selected  to  be  80  pm,  175  pm  and  325  pm 
in  the  present  study.  Although  D16,  D50  and  D84  of  atomized  powders 
strongly  depend  on  the  processing  parameters,  such  as  atomization  gas 
pressure,  superheat  temperature  and  metal/gas  flow  ratio,  the  selection  of  a 
particular  powder  size  distribution  does  not  affect  the  primary  conclusions 
which  will  be  drawn  later.  Taking  into  account  the  maximum  droplet  size 
(~500  pm)  and  the  minimum  droplet  size  (~20  pm)  that  are  generally 
encountered  during  atomization,  the  droplet  sizes  20  pm,  50  pm,  80  pm,  175 
pm,  325  pm  and  500  pm  are  selected  for  calculation  in  the  present  study. 
Following  assumptions  (1)  and  (4)  in  Section  5.2,  representative  initial 
positions  ro=  0  m,  0.00001  m,  0.0001  m,  0.001  m,  0.002  m  and  0.003  m  are 
chosen  for  numerical  calculation  (Zq  remains  to  be  zero  for  all  values  of  ro). 
Therefore,  the  diameter-initial  position  pairs  will  be  {20pm,  (Om,  Om)}, 
{20pm,  (Om,  0.00001m)},  {20pm,  (Om,  O.OOOlm)},  {20pm,  (Om,  0.001m)},  {20pm, 
(Om,  0.002m)},  {20pm,  (Om,  0.003m)},  ...,  {500pm,  (Om,  0.002m)},  {500pm,  (Om, 
0.003m)}. 

No  experimental  values  are  available  for  some  of  the  physical  properties  of 
y-TiAl.  In  this  case,  a  simple  rule  of  mixtures  approach  was  used. 

The  heat  of  fusion  of  y-TiAl  is  determined  on  the  basis  of  information  in 
reference  [109]  as  follows: 
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Y  Ti  A1 

H^J/mol)  H:(J/mol)  H.  (J/mol) 

— ^ - -  0.52— - +  0.48— - +  5.76(J/mol  •  K)  (5-21) 

tL(K)  Ti;(K)  T^(K) 

where  the  superscripts  y,  Ti,  A1  denote  the  properties  of  y-TiAl,  Ti,  and  Al, 
respectively;  and  subscript  m  denotes  melting  temperature. 


The  solid/liquid  interfacial  energy  of  y-TiAl  is  calculated  using  [110]: 

ain(J/m2)=0.46Hfa/atom)[Q(m3/atom)]-2/3  (5-22) 


The  experimentally  determined  and  calculated  physical  properties  of  y-TiAl 
used  in  the  present  study  are  summarized  in  Table  5-1. 


5.4.  Results  and  Discussion 

5.4.1.  Velocity  evolution  of  gas  and  droplets  in  spray  cone 

Figure  5-4  shows  the  axial  gas  velocity  and  droplet  velocities  for  different 

droplet  sizes  along  the  axial  line  (ro=0  m)  as  a  function  of  the  axial  distances. 
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decreasing  droplet  size.  This  is 
consistent  with  the  results 

0 

anticipated  from  Eqns.  (5-5)  and  (5-7),  * 

which  suggest  that,  the  smaller  a  Figure  5-4.  Axial  gas  velocity  and 
droplet  is,  the  larger  the  magnitude  velocities  with  different  droplet 

of  the  force  acted  on  it  becomes.  It  is  sizes  are  20  pm,  50  pm,  80  pm,  175  pm, 
worth  noting  that  the  present 
problem  formulation  assumes  that 

there  is  no  radial  velocity  component  when  ro=0  m,  i.e.,  it  regresses  to  the  one 
dimensional  problem  when  ro=0  m.  Actually,  the  results  summarized  in 
Figure  5-4  are  qualitatively  similar  to  the  one  dimensional  results  obtained  in 
earlier  studies  [88, 92]. 
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Figure  5-5  shows  the  calculated  axial  (Figure  5-5(a))  and  radial  (Figure  5-5(b)) 


components  of  droplet  velocities  and  gas 
velocities  with  ro=0.001  m.  Only  results  for 
20  |im  and  500  pm  droplets  are  shown 
here.  Comparison  of  Figure  5-5(a)  to 
Figure  5-4  reveals  that  both  the  gas  velocity 
and  the  maximum  axial  droplet  velocity 
(for  the  same  droplet  size)  in  Figure  5-5(a) 
are  smaller  than  those  in  Figure  5-4. 
Detailed  calculation  reveals  that,  the 
maximum  axial  component  of  gas  and 
droplet  velocities  decreases  with  increasing 
ro.  Another  distinct  difference  between  the 
results  shown  in  Figure  5-4  and  Figure  5- 
5(a)  is  the  gas  velocities.  It  should  be  noted 
here  that  the  gas  velocities,  both  in  Figure 
5-4  and  Figure  5-5,  are  referred  to  as  the 
velocity  of  gas  in  the  vicinity  of  droplets. 
Accordingly,  although  the  entire  gas 
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Figure  5-5.  Gas  and  droplet 
velocities  for  20  pm  and  500  pm 
droplets  with  an  initial  radial  position, 
ro,  of  0.001  m:  (a)  axial  component 
and  (b)  radial  component. 


velocity  profile  in  the  spray  cone  is 

independent  of  ro,  the  gas  velocities  in  the  vicinity  of  the  droplets  may  be 
different  if  the  trajectories  of  these  droplets  are  different  than  each  other.  It  is 
evident  from  Eqn.  (5-3)  that  the  axial  component  of  gas  velocity  is  a  function 
of  r  and  Z,  i.e.,  the  trajectory  of  droplet.  For  ro=0  m,  the  trajectories  for 


different  droplets  fall  into  one  line,  i.e.,  the  central  line  of  the  spray  cone,  as  a 
result  of  the  absence  of  a  radial  velocity  component  in  this  case.  Therefore, 


when  ro=0  m,  the  axial  component  of  gas  velocity  remains  the  same  for  all 


droplets  (Figure  5-4).  When  ro  deviates  from  zero,  however,  the  trajectories 
of  droplets  with  various  diameters  become  different  from  each  other  (refer  to 
Section  5.4.2),  and  consequently  the  gas  velocities  in  the  vicinity  of  different 


droplets  will  be  different  from  each  other  (Figure  5-5(a)).  As  will  be  shown  in 
Section  5.4.2,  at  any  given  axial  distance,  the  radial  distance  corresponding  to  a 
500  pm  droplet  is  smaller  than  that  of  a  20  pm  droplet.  This  calculated  result, 
along  with  Eqn.  (5-3),  imply  that  the  gas  velocity  in  the  vicinity  of  a  500  pm 
droplet  should  be  larger  than  that  of  a  20  pm  droplet.  The  velocity  profiles  in 
Figure  5-5(b)  are  almost  the  same  as  those  shown  in  Figure  5-5(a).  It  is  worth 
noting,  however,  that  the  radial  velocity  values  for  both  droplet  and  gas  are 
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much  smaller  than  the  axial  ones.  For  example,  the  maximum  radial 
velocity  achieved  by  a  20  pm  droplet  is  4.6  m/s,  while  the  axial  one  is  70  m/ s. 
This  difference  is  qualitatively  consistent  with  available  experimental 
observations  [86],  in  which  the  radial  component  of  droplet  velocity  is  found 
to  be  smaller  than  the  axial  one.  Moreover,  the  radial  component  of  gas 
velocity  in  the  vicinity  of  a  500  pm  droplet,  contrary  to  the  case  of  the  axial 
component,  is  smaller  than  that  in  the  vicinity  of  a  20  pm  droplet  at  any 
given  axial  distance.  This  may  be  understood  from  inspection  of  Eqns.  (5-3) 
and  (5-4),  in  which  'V)rg(r,  Z)  is  a  function  whose  value  increases  with 
increasing  r  provided  that  r  is  not  too  large,  while  Z)  is  a  function 

whose  value  decreases  with  increasing  r.  Regarding  the  influence  of  the 
value  of  ro  on  the  radial  component  of  gas  and  droplet  velocities,  it  is  found 
that,  the  maximum  value  of  the  radial  component  initially  increases  and 
then  decreases  with  increasing  ro. 


5.4.2.  Droplet  size  distribution  in  spray  cone 

Figure  5-6  shows  the  trajectories  of  droplets  of  various  diameters  with 
ro=0.0001  m  and  0.003  m.  The  general  trends  observed  for  ro=0.003  m  may  be 
summarized  as  follows.  First,  for 
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a  given  axial  distance  which  is 

smaller  than  some  critical  value 

(say,  0.75  m),  the  radial  distance 

traveled  by  a  droplet  increases 

with  decreasing  droplet  diameter. 

Second,  for  a  given  axial  distance 

which  is  larger  than  the  critical 

value,  the  radial  location  of 

droplets  becomes  random.  The 

results  obtained  herein  for  all  Figure  5-6.  Trajectories  of  droplet  of  different 
.  .  j .  .  ,  _  sizes  for  ro=0.0001  m  and  ro=0.003  m.  Droplet 

other  initial  conditions  (ro=0,  diameters  are  20  pm,  50  pm,  80  pm,  175  pm, 

0.00001,  0.0001,  0.001  and  0.002  m)  325  pm  and  500  pm. 
are  qualitatively  similar  to  those 

obtained  with  ro=0.003  m.  It  is  worth  noting,  however,  that  the  radial 
distance  range,  Ar,  at  a  given  axial  distance  decreases  with  decreasing  ro.  For 
example,  when  ro=0.0001  m,  the  radial  distance  range  Ar  corresponding  to  any 
axial  distance  is  so  small  that  the  trajectories  of  different  droplets  are  almost 
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indistinguishable  in  Figure  5-6.  In  addition,  for  the  same  droplet  size, 
droplets  initially  located  at  a  larger  ro  will  go  farther  in  radial  direction  from 
the  spray  cone  axis  than  that  at  a  smaller  ro- 

For  an  individual  droplet,  the  radial  distance  corresponding  to  a  given  axial 
distance  may  be  expressed  as  an  integral  of  the  product  of  flight  time  and  the 
radial  droplet  velocity  component.  Consequently,  there  will  be  several  factors 
influencing  the  radial  distances  traveled  by  droplets  at  a  given  axial  distance, 
which  finally  leads  to  the  calculated  two  dimensional  droplet  size  distribution 
in  the  spray  cone.  For  droplets  with  identical  ro,  the  radial  droplet  velocity 
component  for  a  fine  droplet  is  initially  larger  and  then  smaller  than  that 
corresponding  to  a  coarse  one  (Figure  5-5(b)).  Moreover,  as  will  be  discussed 
in  a  later  section,  the  flight  time  for  a  given  axial  distance  critically  depends 
on  droplet  size.  Combined,  all  of  these  effects  result  in  the  trajectory  band  for 
each  ro,  as  shown  in  Figure  5-6.  The  effect  of  ro  on  the  radial  distance  traveled 
by  a  droplet  may  be  attributed  to  two  factors.  First,  to  reach  a  given  axial 
distance,  droplets  which  are  initially  located  at  a  larger  ro  will  experience  a 
longer  flight  time  relative  to  those  corresponding  to  a  smaller  ro,  as  will  be 
shown  in  Section  5.4.3.  Second,  as  discussed  in  Section  5.4.1,  the  radial 
component  of  droplet  velocity  initially  increases  and  then  decreases  with 
increasing  ro.  Detailed  calculations  for  other  initial  positions  (ro=0,  0.00001, 
0.001  and  0.002  m)  indicate  that  the  radial  distance  traveled  by  a  droplet  is 
more  significantly  affected  by  changes  in  flight  time  than  that  in  radial 
component  of  droplet  velocity. 

It  may  be  envisioned  that,  if  all  the  results  calculated  for  different  values  of  ro 
were  to  be  shown  in  Figure  5-6,  the  trajectory  bands  for  different  values  of  ro 
will  overlap  with  each  other  and  a  complete  picture  of  the  spatial  distribution 
of  droplets  with  different  diameters  in  the  spray  cone  will  be  obtained.  As  a 
result,  when  the  axial  distance  is  smaller  than  some  critical  value  (say  0.75  m), 
coarse  droplets  tend  to  concentrate  in  the  central  region  of  the  spray  cone, 
whereas  fine  droplets  populate  the  periphery;  when  the  axial  distance  is  larger 
than  the  critical  value,  the  two  dimensional  droplet  size  distribution  in  the 
spray  cone  will  be  homogeneous.  These  calculated  results  are  qualitatively 
consistent  with  the  experimental  results  reported  by  McDonell  et  al.  [111]. 
The  experimental  results  of  McDonell  et  al.  [Ill]  reveal  that  coarse  droplets 


concentrate  in  the  central  region  of  the  spray  cone  at  an  axial  distance  of  0.025 
m.  When  the  axial  distance  increases  from  0.025  m  to  0.28  m,  however,  the 
two  dimensional  droplet  size  distribution  becomes  relatively  insensitive  to 
the  radial  position,  i.e.,  from  heterogeneous  to  almost  homogeneous.  The 
quantitative  discrepancy  between  the  results  obtained  herein  and  those 
obtained  by  McDonell  et  al.  [Ill]  may  be  attributed  to  the  difference  in 
material  properties  and  process  parameters.  The  study  by  McDonell  et  al. 
[Ill],  for  example,  involved  the  atomization  of  water  with  air. 

5.4.3.  Flight  time 

The  flight  times  corresponding  to  droplets  with  different  sizes  and  initial 
conditions  were  calculated  as  a  function  of  the  axial  distance.  The  calculated 
results  reveal  that,  for  the  same  values  of  ro,  there  are  two  distinct  regimes. 
For  axial  distances  smaller  than  a  critical  value  (e.g.,  0.35  m  for  ro=0.0001  m), 
flight  time  increases  with  increasing  droplet  size.  This  may  be  rationalized  on 
the  basis  of  the  results  shown  in  Figures  4  and  5(a),  which  reveal  that  fine 
droplets  are  readily  accelerated  and  that  the  axial  velocity  component  of  fine 
droplets  exceeds  that  of  coarse  droplets.  Accordingly,  the  time  elapsed  during 
travel  of  a  particular  distance  decreases  with  decreasing  droplet  size.  When 
the  axial  distance  is  larger  than  the  critical  value,  the  flight  time 
corresponding  to  finer  droplets  (i.e.,  20  (tm,  50  |xm  and  80  |im)  increases 
rapidly  to  a  much  higher  magnitude  relative  to  that  corresponding  to  coarse 
droplets  (i.e.,  175  |im,  325  pm  and  500  pm).  This  behavior  may  be  attributed 
to  the  fact  that  finer  droplets  are  readily  decelerated  to  very  small  velocity 
values  once  a  maximum  velocity  has  been  achieved,  while  the  velocities  of 
coarse  droplets  remain  almost  unchanged  once  the  maximum  velocity  has 
been  attained  (e.g.,  see  Figures  4  and  5).  When  ro  increases,  flight  time 
increases  rapidly  due  to  the  differences  in  the  magnitude  of  the  axial  velocity 
component  of  the  droplets,  as  discussed  in  detail  in  Section  5.4.1. 

5.4.4.  Cooling  rate 

The  cooling  rates  of  a  droplet  in  various  thermal  and  solidification  stages  are 
completely  different  than  each  other,  from  being  negative  in  recalescence 
stage  to  being  zero  in  peritectic  solidification  stage,  and  to  being  positive  in 
other  three  stages  (cooling  in  the  liquid  state,  cooling  in  the  solid  state  and 
segregated  solidification).  Even  in  the  same  stage,  cooling  rate  of  a  droplet 
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also  critically  depends  on  its  droplet  size,  initial  position  and  axial  position. 
For  example,  the  cooling  rate  in  the  stage  of  cooling  in  the  liquid  state  ranges 
from  1.73x103  K/s  to  1.55  xlO^  K/s,  depending  on  its  droplet  size,  initial 
position  and  axial  distance  from  the  atomization  zone.  Detailed  calculation 
of  cooling  rates  for  droplets  of  different  sizes  and  different  initial  conditions, 
however,  indicates  that  there  are  two  distinct  tendencies:  i)  cooling  rate 
decreases  significantly  with  increasing  droplet  size;  and  ii)  cooling  rate 
remains  almost  unchanged  when  ro  increases.  As  evident  from  Eqn.  (5-10), 
the  difference  between  the  calculated  cooling  rate  for  different  values  of  D 
results  from  the  variations  in  droplet  volume,  V,  and  heat  transfer 
coefficient,  h.  The  latter,  h,  is  a  complex  function  of  D  and  Udg,  the  relative 
velocity  between  gas  and  droplet.  The  influence  of  the  value  of  ro  on  the 
calculated  cooling  rate  may  be  inferred  from  the  variation  of  Udg  under 
different  ro,  as  discussed  in  Section  5.4.1. 
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Figure  5-7.  Thermal  history  of  different  droplets  with  initial  radial  position 
ro=0  m;  (a)  D=20  pm;  (b)  D=80  pm;  (c)  D=175  pm  and  (d)  D=325  pm. 

5.4.5.  Thermal  history 

Figure  5-7  shows  the  calculated  thermal  histories  for  droplets  of  different  sizes 
assuming  a  value  of  ro=0  m.  All  five  stages  described  in  Section  5.2.3  were 


experienced  by  droplets  with  D=20  pm  (Figure  5-7(a))  and  D=  80  pm  (Figure  5- 
7(b)).  Some  stages,  such  as  recalescence  and  cooling  in  the  solid  state,  are 
absent,  depending  primarily  on  droplet  size.  The  disappearance  of  a 
recalescence  stage  in  coarse  droplets  (D=175  pm  and  325  pm)  may  be 
rationalized  on  the  basis  of  the  very  small  undercoolings  that  are  associated 
with  these  droplets  (see  Eqn.  (5-9)).  The  disappearance  of  the  cooling  in  solid 
state  stage,  as  noted  for  325  pm  droplets,  is  a  result  of  their  relatively  lower 
cooling  rate  compared  with  those  of 
finer  droplets  (see  Section  5.4.4). 

Accordingly,  in  large  droplets,  the 
temperature  is  not  sufficiently  low  in 
the  calculated  regime  to  develop  this 
stage.  Detailed  calculation  indicates 
that,  when  ro  increases,  the  curve 
becomes  steeper.  This  phenomenon 
is  a  result  of  the  combined  effects  of 
the  longer  flight  time  and  the 
relatively  unchanged  cooling  rate 
with  increasing  ro- 
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5.4.6.  Solidification  behavior 
Figure  5-8  shows  the  solidification 
tendencies  of  droplets  of  different 

sizes  for  two  different  initial 

A  n  nnnni  .v,  ,^..4  n  no's  Figure  5-8.  Solidification  tendencies  of 

conditions  (ro=0.00001  m  and  0.003  ,^,^0  00001  m 

m).  The  results  obtained  for  other  and  (b)  ro=0.003  m.  The  droplet  diameters 

initial  conditions  (ro=0  m,  0.0001  m,  ^  ^0  50  mn,  80  urn,  175  nm,  325 

0.001  m  and  0.002  m)  are  similar  to 

those  shown  in  Figure  5-8.  There  are  two  trends  that  are  evident  from  the 
results  summarized  in  Figure  5-8.  The  first  one  is  that  for  some  axial  distance 
(say  0.14  m  in  Figure  5-8(a)),  the  fraction  solidified  in  coarse  droplets  (0.13  in 
175  pm  droplets)  is  higher  than  that  in  fine  droplets  (0  in  80  pm  droplets). 
The  second  one  is  that,  once  the  nucleation  temperature  is  reached  in  fine 
droplets  (80  pm),  they  are  more  readily  solidified  relative  to  coarse  droplets 
(175  pm)  even  if  they  start  to  nucleate  much  later  than  coarse  ones.  The  first 
trend  may  be  rationalized  from  the  very  small  actual  undercoolings 
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experienced  by  coarse  droplets,  which  result  in  a  much  higher  nucleation 
temperature  relative  to  those  experienced  by  fine  droplets.  Consequently, 
although  fine  droplets  have  higher  cooling  rate  than  coarse  ones,  it  is  still 
possible  that  in  some  cases  the  nucleation  temperature  of  coarse  droplets  is 
reached  earlier  than  that  of  fine  droplets,  i.e.,  coarse  droplets  start  nucleation 
earlier  than  fine  ones.  The  second  trend  may  be  understood  from  inspection 
of  Eqns.  (5-11),  (5-14),  (5-19b)  and  (5-20).  For  example,  Eqn.  (5-14),  which 
describes  the  solidification  occurring  in  the  recalescence  stage,  implies  that  a 
smaller  droplet  size  and  a  larger  undercooling  correspond  to  a  higher 
solidification  rate.  A  smaller  droplet  size  also  corresponds  to  a  larger 
undercooling  (Table  5-2).  Therefore,  fine  droplets  solidify  at  a  higher  rate 
relative  to  coarse  ones  during  recalescence  stage.  Furthermore,  the  second 
trend  is  also  consistent  with  the  results  obtained  from  the  one  dimensional 
modeling  of  spray  atomization  [88]. 

Comparison  between  Figure  5-8(a)  and  Figure  5-8(b)  indicates  that,  for 
identical  droplet  sizes,  the  solidification  fraction  at  a  given  axial  distance 
increases  with  increasing  ro.  This  may  be  understood  in  terms  of  the  flight 
time  which  increases  rapidly  with  increasing  ro,  as  discussed  in  Section  5.4.3. 
A  longer  flight  time  will  allow  for  more  fraction  being  solidified  provided 
that  the  droplets  under  investigation  are  not  fully  solidified. 

The  two  dimensional  solidification  behavior  of  atomized  y-TiAl  droplets  in 
the  spray  cone  may  be  elucidated  through  calculations  of  the  average  solid 
fraction,  favg/  for  a  collection  of  droplets  which  is  initially  located  at  an 
identical  position,  ro.  In  a  given  collection  of  droplets,  the  average  solid 

fraction  in  it  may  be  estimated  from  [27]: 

favg  =  (fl6  +  f84)  X  32%  +  fso  X  36%  (5-23) 

where  fi6,  fso  and  f84  represent  the  solid  fraction  in  characteristic  droplets  with 

diameters  Di6,  D50  and  D84,  respectively. 

The  size  distribution  of  a  collection  of  droplets  with  identical  initial  position, 
ro,  may  be  assumed  to  be  the  same  as  that  corresponding  to  the  droplets 
formed  in  the  atomization  zone,  provided  that  droplets  are  randomly 
distributed  in  the  atomization  zone  once  they  are  formed  (refer  to 
assumption  (4)  in  Section  5.2).  Accordingly,  Die,  D50  and  D84  in  a  collection  of 


droplets  with  a  particular  ro  will  be  the  same  for  all  collections  of  droplets 
with  different  given  initial  positions,  and  representative  of  the  size 
distribution  of  droplets  formed  in  the  atomization  zone,  which  are  80  |xm,  175 
pm  and  325  pm,  respectively.  Hence,  favg  for  a  collection  of  droplets  with  any 
particular  ro  may  be  obtained  by  using  Eqn.  (5-23)  along  with  the  magnitude  of 
fi6,  fso  and  f84  in  droplets  with  diameters  of  80  pm,  175  pm  and  325  pm, 
respectively.  It  should  be  noted  that  the  solid  fractions  fi6,  fso  and  f84  are 
functions  of  the  initial  position,  ro. 
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The  calculated  favg  as  a  function  of 

axial  distance  is  summarized  in 

Figure  5-9  for  ro=0,  0.001  and  0.003 

m.  The  effect  of  changing  ro  on  the 

solidification  behavior  is  evident 

there:  the  larger  the  value  of  ro,  the 

higher  the  solid  fraction  at  a  given 

axial  distance.  As  mentioned  in 

Section  5.4.2,  for  identical  droplet 

size  (D),  the  radial  distance  traveled 
u  j  Figurc  5-9.  The  average  solid  fraction 

by  a  drop  y  g  contained  in  droplets  which  are  initially  located 

distance  increases  with  increasing  at:  (a)  ro=0  m;  (b)  ro=0.001  m  and  (c) 

ro,  i.e.,  if  a  substrate  is  inserted  ro-0.003  m. 

across  the  spray  cone,  droplets  with  smaller  ro  will  concentrate  in  the  central 
region  of  the  substrate  whereas  droplets  with  larger  tq  will  populate  the 
periphery  of  the  substrate.  Therefore,  if  a  substrate  were  to  be  inserted  across 
the  spray  cone,  the  solid  fraction  contained  in  the  droplets  arriving  on  the 
substrate  would  become  higher  with  increasing  radial  distance  from  the 
center  of  the  substrate.  Although  there  is  no  direct  experimental  support  to 
the  two  dimensional  solidification  behavior  calculated  herein,  it  is 
experimentally  found  that  in  spray  deposited  materials,  such  as:  Al-Cu-Li-Ag- 
Mg-Zr  alloy  WeldaliteTM  049  [112],  6061  aluminum  alloy  [113]  and  more 
recently  for  spray  formed  y-TiAl  [114],  the  peripheral  region  generally 
contains  a  higher  volume  fraction  of  porosity  relative  to  that  present  in  the 
central  region.  Moreover,  the  available  experimental  evidence  suggests  that 
porosity  present  in  spray  deposited  material  generally  develops  from 
interstices  formed  as  droplets  impinge  on  one  another,  leaving  micro-sized. 


irregular  cavities  as  they  overlap  [112,  115].  On  the  basis  of  this  mechanism  of 
porosity  formation  in  spray  deposited  materials,  a  high  solid  fraction  in  the 
spray  will  result  in  a  high  volume  fraction  of  porosity  in  the  deposit,  since 
insufficient  molten  material  is  available  to  flow  over  the  deposit  to  fill  pores 
and  interstices.  Therefore,  the  experimental  observation  that  the  peripheral 
region  contains  a  higher  volume  fraction  of  porosity  relative  to  that  present 
in  the  central  region  implies,  to  some  extent,  that  the  fraction  solidified  in  the 
periphery  is  higher  than  that  in  the  central  region. 

5.5.  Summary 

1.  The  velocity,  temperature,  cooling  rate,  flight  time  and  solidification 
behavior  of  atomized  droplets  strongly  depends  on  the  initial  position,  ro,  and 
the  diameter  of  the  droplet,  D.  In  addition,  the  velocity  profiles,  temperature 
histories,  cooling  rates,  flight  times  and  solidification  behavior  of  droplets  are 
also  closely  related  to  each  other. 

2.  For  identical  ro,  i)  fine  droplets  are  more  readily  accelerated  or  decelerated 
relative  to  coarse  droplets  along  both  axial  and  radial  directions,  ii)  fine 
droplets  cool  at  a  higher  cooling  rate  than  coarse  ones,  and  iii)  fine  droplets 
also  solidify  at  a  faster  rate  relative  to  coarse  ones. 

3.  For  identical  D,  i)  the  maximum  value  of  the  axial  component  of  droplet 
velocity  decreases  with  increasing  ro,  while  the  maximum  value  of  the  radial 
component  of  droplet  velocity  initially  increases  and  then  decreases  with 
increasing  ro,  ii)  flight  time  increases  with  increasing  ro  and  iii)  the  solid 
fraction  at  any  given  axial  distance  increases  with  increasing  ro  resulting  from 
the  longer  flight  time  for  larger  ro. 

4.  The  two  dimensional  droplet  size  distribution  in  the  spray  cone  changes 
from  being  heterogeneous  to  be  almost  homogeneous  with  increasing  axial 
flight  distance.  When  the  axial  distance  is  smaller  than  a  critical  value 
(approximately  0.75  m  for  the  conditions  used  in  present  study),  coarse 
droplets  concentrate  in  the  central  region  of  the  spray  cone,  whereas  fine 
droplets  populate  the  periphery.  When  the  axial  distance  increases  beyond 


this  critical  value,  the  spatial  droplet  size  distribution  in  the  spray  cone 
becomes  independent  of  the  radial  position. 

5.  The  two  dimensional  distribution  of  the  fraction  solidified  in  the  spray 
cone  is  heterogeneous.  The  fraction  solidified  at  any  axial  distance  increases 
with  increasing  radial  distance  from  the  spray  axis. 


CHAPTER  6. 


FINITE  ELEMENT  SIMULATION  OF  CREEP  DEFORMATION 
AND  RUPTURE  OF  TITANIUM  ALUMINIDES 


6.1.  Objective 

The  objectives  of  this  work  are:  1)  to  study  the  creep  behavior  of  TiAl  alloys 
under  tensile  creep  conditions  and  to  investigate  the  effect  of  grain  boundary 
sliding  on  the  overall  high  temperature  deformation  and  damage  processes, 
and  to  compare  the  creep  behavior  of  FL  and  duplex  microstructures;  2)  to 
model  and  analyze  the  grain  boundary  cavitation  in  fully  lamellar  (a2+y) 
titanium  aluminide. 


6.2.  Model  Formulation 


6.2.1.  Creep  deformation 

The  FL  and  duplex  forms  of  TiAl  are  modeled  using  a  nonlinear,  implicit  two 

dimensional  finite  element  code,  NIKE2D  [116],  developed  at  Lawrence 

Livermore  National  Laboratory. 

The  representative  microstructure 

is  taken  as  a  conglomeration  of 

seven  equiaxed  grains.  The  duplex 

structure  is  modeled  by 

considering  a  single  equiaxed  grain 

of  a2,  surrounded  by  six  equiaxed  y 

grains.  Grains  of  the  FL 

microstructure  consist  of  alternate 

a2  laths  and  packets  of  y  laths  as 

shown  in  the  Figure  6-1.  The 

orientation  of  the  laths  with 

respect  to  each  other  are  chosen  at 

random  so  that  the  model 

resembles  a  realistic  lamellar 

TiAl  microstructure.  The  darker  shading 

microstructure  [117, 118].  indicates  the  a2  lamellae. 


Two  case  studies  have  been  performed  for  both  the  duplex  and  FL  models.  In 
the  first  case,  the  lath  and  grain  boundaries  are  not  allowed  to  slide.  In  the 
second  case,  all  grain  boundaries,  except  those  at  the  continuum  boundaries, 
for  both  models  are  capable  of  frictionless  sliding.  The  results  of  the  present 
model  have  been  compared  with  those  of  a  model  having  fifteen  grains  and  it 
has  been  found  that  the  stress  distributions  are  identical.  This  finding 
provides  justification  for  using  the  present  smaller  model  which  reduces  the 
computation  time  without  adversely  affecting  the  accuracy  of  the  analysis. 
Four  noded  isoparametric  elements 
are  used  to  model  the  duplex  and  FL 
forms.  Sliding  elements  are  used 
for  modeling  the  grain  boundaries. 

A  penalty  fimction  approach  is  used 
in  NIKE2D  to  implement  the 
contact  mechanics  formulation  that 
takes  into  account  of  the  sliding 
capability  of  the  elements.  A  2x2 
gauss  quadrature  is  used  for 
numerical  integration  and  the 
Broyden-Fletcher-Goldfarb-Shanno 
(BFGS)  method  [119]  is  employed  as 
the  solution  routine  for  the  set  of 
non-linear  equations  arising  out  of 
the  finite  element  formulation. 

6.2. 1.1.  Geometry  of  the  model 
The  geometry  of  the  model  is  shown  in  Figure  6-2.  Both  the  FL  and  duplex 
forms  are  assumed  to  have  a  grain  size  of  70  mm  and  have  a  phase  ratio  of  1:3 
(a2 :  Y)-  The  orientation  of  the  grains  for  FL  models  are  given  in  Tables  6-1 
and  6-2: 

6.2. 1.2.  Boundary  conditions 

The  boundary  conditions  of  the  present  models,  illustrated  in  Figure  6-2,  are 
implemented  in  such  a  way  that  the  grain  boundaries  can  slide  freely  without 
separation.  Nodes  at  the  grain  boundary  triple  points  are  pinned  together  to 
meet  strain  compatibility  conditions.  Symmetry  boundary  conditions,  u=0 


2 

Figure  6-2.  Unit  cell  with  grain  boundary 
conditions. 


and  v=0,  are  applied  on  sides  OF  and  OA  of  the  cell  OACF  respectively.  This 
unit  cell  forms  a  larger  cell  when  reflected  about  OA  and  OF  and  repetition  of 
this  process  builds  a  representative  bulk  material  sample.  Nodes  on  cell 
boundary  AC  are  constrained  to  have  the  same  y  displacement,  u,  while 
nodes  on  boundary  CF  are  constrained  to  have  the  same  z  displacement,  v,  at 
all  times.  Stress,  E,  is  applied  on  boundary  CF.  External  stress  is  maintained 
constant  all  the  time.  Plane  strain  conditions  are  assumed  to  exist  in  the 
representative  continuum  [117,  118]. 


Table  6-1.  Angular  lath 
orientation  in  the  grains 
of  the  FL-1  model. 


Grain 

no. 

Orientatio 
n,  (p 

(degrees) 

1 

10.0 

2 

125.2 

3 

99.1 

4 

49.3 

5 

99.1 

6 

10.0 

7 

125.2 

Table  6-2.  Angular 
lath  orientation  in  the 
grains  of  the  FL-2 
model. 


Grain 

Orientatio 

no. 

n,  (p 

(degrees) 

1 

10.0 

2 

125.2 

3 

99.1 

4 

49.3 

5 

19.6 

6 

169.5 

7 

8.8 

6.2.2.  Creep  constrained  cavitation 

At  very  low  tensile  stresses,  grain  boundary  diffusion  is  the  dominant 
mechanism  and  creep  flow  is  negligible  [120].  The  volumetric  growth  rate  of 
equilibrium  shaped  cavities  in  this  case  is  given  after  Rice  and  Needleman 
[120]  as: 


(Tn-(l-C)(ys 

V,  =47t(D - 

'  log(l/0-(3-0(l-0/2 


(6-1) 


Where,  On  is  the  average  normal  stress  on  the  cavitated  grain  boundary.  The 
sintering  stress,  Os,  is  given  as  as=2  Ys  sin(\|r)/a,  where  js  is  the  surface  free 
energy.  The  cavity  tip  angle,  \|/,  is  obtained  from  the  following  relation: 

Cos(i//)  =  ^  (6-2) 


where  yb  is  the  grain  boundary  free  energy.  The  area  fraction  of  the  cavitated 


grain  boundary,  is  given  by:  . 

^  =  maximum^ 
where  b  is  the  cavity  half  spacing 


a 


f-). 

yb)  ya  +  1.5L,  _ 
as  shown  m  Figt 


re  6-3.  L  is  given  by: 


(6-3) 

Here,  ae  is  the  remote  effective  Mises  stress  and  is  the  corresponding 
effective  creep  strain  rate  which  are  obtained  from  the  power  law. 


D  is  defined  by: 


‘D  =  Dg8^aikT 


where  Db  is  the  grain  boundary  diffusivity,  8b  is 
the  grain  boundary  thickness.  Omega  is  the  atomic 
volume  and  kT  is  the  energy  per  atom  measure  of 
temperature.  L  denotes  the  effective  radius  of  the 
grain  boundary  diffusion  zone  around  the  cavity. 
Material  outside  this  zone  primarily  deforms  solely 
by  dislocation  creep  while  within  this  zone  grain 
boundary  diffusion  from  the  cavity  surface 
contributes  to  the  deformation. 


(6-4) 

Applied  Stress  Direction 

k 


2Ro 


2a 


2b, 


Applied  Stress  Direction 

Figure  6-3.  Schematic 
showing  grain  boundary 
cavity. 


Equations  for  void  growth  due  to  dislocation  creep 
were  developed  by  Budiansky  et  al.  [121]  for  high 

triaxial  stress  conditions.  Tvergaard  [122]  modified  these  equations  to  make 
them  applicable  for  the  low  triaxiality  range  and  have  the  form: 


V^=+2Ke„a%{\if) 


and 


a. 


+  Pn 


for-i-^>l 

-  <y. 


V,  =  a^h{xi/)[a„+P„]^  for 


(6-5) 


<1 


where  an=3/2n  and  Pn=(n-l)(n+0.4319)/(n2).  The  parameters  Om  and  Og  are 
the  remote  mean  stress  and  Mises  stress  respectively.  The  function  h(\}/) 
accounts  for  the  shape  of  the  cavity  which  is  given  as: 

h(  y/)  =  [(1  +  cos(  v/-))”'  -  iCos(  v^)]  /  Sin(  iff)  (6-6) 

The  above  equation  is  based  on  the  assumption  that  surface  diffusion  in  the 
voids  is  sufficiently  rapid  to  retain  a  spherical  cap  (equilibrium)  shape. 
Cavities  grow  as  a  result  of  diffusion  of  surface  atoms  of  the  cavity  to  adjacent 
grain  boundaries.  If  the  void  growth  is  facilitated  by  creep  deformation,  the 
cavities  will  tend  to  elongate  in  the  tensile  direction.  On  the  other  hand,  the 


cavities  tend  to  become  lenticular  when  grain  boundary  diffusion  is  the  only 
mechanism  of  void  growth.  It  is  rapid  surface  diffusion  combined  with  grain 
boundary  diffusion  and  dislocation  creep  that  results  in  equilibrium  shaped 
cavities  [120].  The  total  volumetric  cavity  growth  rate  can  be  determined  by 
summing  equations  (6-1)  and  (6-5),  thus: 

V  =  V,-tV2,for— <10  (6-7) 

L 

The  corresponding  rate  of  change  of  cavity  radius  is  given  as: 

a  = - -  (6-8) 

Am  h{  yr) 

The  growth  of  cavities  on  grain  boundary  facets  results  in  the  separation  of 
grains  by  the  plating  of  atoms  out  onto  the  grain  boundary  from  the  surface  of 
the  cavities  [123]. 


The  average  separation  delta  can  be  expressed  in  terms  of  cavity  volume  V 
and  cavity  half  spacing  b  as  5  =  V/irb^.  Therefore,  the  rate  of  separation  of  the 
grains,  6,  is  [122]: 


3  = 


V  2Vb 


(6-9) 


Ttb^  Ttb^ 

where  V  is  given  by  equation  (6-7).  The  contribution  of  the  second  term  in 
equation  (6-9)  becomes  important  when  finite  strains  in  the  zones  close  to  the 
cavities  result  in  significant  changes  in  the  void  spacing.  The  rate  of  change 
of  void  spacing,  b,  can  be  correlated  with  the  true  strain  rate,  Sp,  in  the  grain 
boundary  plane  under  plane  strain  conditions  as  in  [124]: 


b 

b 


(6-10) 


We  assume  that  all  the  cavities  along  a  facet  are  of  same  size  and  they  are 
spaced  equally.  In  this  case,  the  distinct  distribution  of  cavities  at  grain 
boundaries  can  be  represented  by  a  uniform  value  of  6  [124]. 


6.2.2. 1.  Method  of  analysis 

It  is  assumed  that  the  spatial  distribution  of  the  cavitating  facets  follows  a 
periodic  order  so  that  we  can  confine  our  analysis  to  a  unit  cell.  Only 
transverse  grain  boundaries  with  respect  to  the  direction  of  externally  applied 
stresses  are  assumed  to  cavitate.  Figure  6-2  shows  one  quarter  of  the  unit  cell 
with  boundary  conditions.  It  can  be  seen  from  this  figure  that  the  dimensions 


of  the  unit  cell  are  specified  in  terms  of  the  number  of  grains  in  the  y 
direction,  mi,  the  number  of  grains  in  the  z  direction  m2  and  the  initial 
radius  of  the  grain  boundary  facets  Rq. 

As  in  the  creep  deformation  models,  finite  element  techniques  were  used 
here  also  to  model  creep  constrained  cavitation  in  the  fully  lamellar  TiAl. 
All  the  analyses  are  performed  under  plane  strain  boundary  conditions. 
Finite  element  code,  NIKE2D  [116]  and  its  associated  pre-  and  post-processors 
were  used  to  analyze  all  of  the  cases  presented  here.  Grain  boundary  sliding 
under  the  action  of  externally  applied  stresses  was  simulated  using  slide  lines 
along  the  grain  boundaries.  Details  on  the  use  of  slide  lines  in  FE  models  can 
be  found  elsewhere  [116,  125].  In  all  of  the  analyses  performed  here,  it  is 
assumed  that  grain  boundary  sliding  takes  place  freely.  This  assumption 
should  be  particularly  valid  for  intermetallic  materials,  such  as  TiAl  alloys, 
which  have  a  relatively  low  dislocation  mobility  under  the  conditions 
considered  making  boundary  sliding  more  prevalent  as  a  deformation 
mechanism  [126,  127].  Furthermore,  Hayes  and  Martin  [9]  have  observed 
grain  boundary  triple  point  wedge  cracks  in  single  phase  gamma  TiAl  which 
suggest  that  the  friction  stress  against  grain  boundary  sliding  in  this  material 
is  relatively  low  under  creep  conditions. 

It  has  also  been  proposed  that  certain  y/y  lath  boundaries  may  slide  by  the 
dissociation  of  superdislocations  in  fully  lamellar  TiAl  [128].  Although 
experimental  evidence  has  been  reported  for  this  mechanism  during  plastic 
deformation  at  ambient  temperature,  both  theoretical  [61]  and  experimental 
[43]  studies  have  shown  that  with  the  possible  exception  of  the  earliest  stages 
of  primary  creep,  this  mechanism  is  not  operative  under  typical  creep 
conditions.  Accordingly,  lath  boundary  sliding  is  not  considered  in  the 
present  models. 

The  boundary  conditions  at  the  cavitating  grain  boundary  are  modeled  after 
Tvergaard  [122].  In  each  time  increment,  the  facet  OH  was  subjected  to  an 
incremental  displacement  5  given  as: 

S 


where  5  is  given  by  equation  (6-11)  and  5t  is  the  time  increment.  Stresses  am, 
<Je  and  an  used  in  equations  (6-1)  and  (6-5)  and  the  true  strain  rate  8p  used  in 
evaluating  equation  (6-10)  are  obtained  directly  from  the  elements  at  the 
center  of  the  cavitating  facet  at  each  time  step.  The  center  of  the  cavitating 
facet  is  chosen  as  this  sampling  point  because  the  influence  of  the  singularity 
at  the  grain  boundary  triple  points  is  minimal  there  and  thus  the  stresses  and 
strains  are  more  representative  of  the  overall  stress  on  the  boundary.  For 
fully  lamellar  models,  this  position  is  also  within  a  y  lath  which  is  under  a 
higher  stress  compared  to  neighboring  a2  laths  since  a2  is  weaker  than  y 
under  the  present  creep  conditions  [61].  The  length  of  the  time  steps  during 
the  initial  period  of  cavitation  is  kept  very  small  so  that  numerical  stability  is 
maintained.  Normal  stresses,  Gn,  at  the  cavitating  grain  boundary  elements 
are  evaluated  at  each  time  step  to  assure  that  they  are  greater  than  zero. 
Whenever,  the  normal  stress  reaches  zero  in  any  element  at  a  cavitating 
boundary,  that  particular  element  edge  is  treated  as  a  free  surface  from  that 
time  step  onward. 

62.2.2.  Modeling  TiAl  microstructures 

Models  of  the  dual  phase  equiaxed  microstructure  that  has  the  same  phase 
ratio  (a2/y)  as  that  in  the  fully  lamellar  microstructure  have  been  developed 
to  compare  the  latter  structure  with  the  former. 

Two  forms  of  dual  phase  a2+y  microstructures  are  analyzed  here.  The  first, 
DX-1,  is  modeled  by  considering  a  unit  cell  which  has  equiaxed  grains  of  a2 
phase  surrounded  by  six  equiaxed  grains  of  gamma  phase.  Thus  the  mesh 
used  here  is  the  same  as  that  used  in  the  single  phase  model  except  that  the 
central  grain,  4,  is  composed  of  a2  phase.  The  boundary  conditions  are 
identical  to  those  in  the  previous  case.  The  second  model  of  the  dual  phase 
structure,  DX-2,  is  configured  as  grains  3  and  5  of  a2  phase  and  all  other 
grains,  including  grain  4,  of  y  phase.  Other  than  this  difference,  the  mesh 
used  for  the  DX-2  model  is  the  same  as  that  used  for  the  DX-1  model. 

Furthermore,  it  has  been  proposed  that  serrated  grain  boundaries  inhibit 
grain  boundary  sliding  during  creep  deformation  [48,  61,  129,  130]. 
Considering  this  hypothesis,  models  of  fully  lamellar  TiAl  intermetallics  are 
analyzed  in  the  present  work  for  both  the  case  of  cavitation  with  the  presence 


of  grain  boundary  sliding  and  the  case  of  cavitation  without  the  influence  of 
grain  boundary  sliding.  We  refer  to  the  mesh,  used  in  creep  deformation 
studies  as  either  FL-1  or  FL-2  depending  on  the  loading  and  boundary 
conditions.  For  the  FL-1  models,  tensile  stress  X  is  applied  at  the  face  CF  of  the 
unit  cell  OACF  and  all  other  boundary  conditions  are  the  same  as  for  the 
single  phase  and  dual  phase  models.  For  the  FL-2  models,  tensile  stress  X  is 
applied  at  the  face  OA  and  the  boundaries  OF  and  CF  are  constrained  to  have 
zero  displacements  along  y  and  z  directions  respectively.  The  nodes  at 
boundaries  OA  are  constrained  to  move  together  along  the  z  direction  and 
the  nodes  at  boundaries  AC  are  constrained  to  move  together  along  the  y 
direction  at  all  time.  As  before,  grain  boundary  triple  points  are  pinned 
together  to  have  the  same  y  and  z  displacements  at  all  the  time.  The 
cavitating  facet  is  OH  in  all  models  except  for  the  FL-2  models  in  which  it  is 
the  facet  EF. 

The  last  stages  of  cavity  coalescence  involves  rapid  tearing  of  the  ligaments 
between  the  cavities  and  it  is  reasonable  to  assume  that  this  process  begins  to 
occur  once  the  cavitation  damage,  co,  reaches  about  0.6  [131].  Accordingly,  a 
critical  level  of  cavitation  is  assumed  to  occur  when  a)=0.6  in  all  of  the  present 
analyses.  In  all  of  the  analyses,  it  is  assumed  that  the  cavities  nucleate  upon 
loading  and  there  is  no  nucleation  of  any  new  cavities  during  the  rest  of  the 
simulation.  The  following  initial  values  are  assumed,  (a/b)o=0.1  and 
(bo/Ro)=0-4.  It  should  be  noted  that  the  above  assumptions  are  uniformly 
made  for  all  of  the  present  models.  Riedel  [131]  observed  that  the  initial  size 
of  the  cavities  in  metals  and  alloys  is  typically  0.2-1  pm.  Accordingly,  we 
assume  a  representative  initial  value  of  a  to  be  0.7  pm.  We  note  that  the 
model  results  are  relatively  insensitive  to  the  initial  cavity  radius  as  long  as  it 
is  within  a  range  of  reasonable  values.  For  most  metals  and  alloys,  the  grain 
boundary  surface  energy,  js  is  between  1-2  J/m^  [132].  We  estimate  ys=pb\10 
=1.69  J/m2,  where  b  is  the  burgers  vector  and  p  is  the  shear  modulus, 
following  the  work  of  Chu  and  Thompson  [132].  Sintering  stress,  Os,  is  found 
to  be  4.54  MPa  for  a  cavity  radius  of  a  =  0.7  pm  which  is  negligibly  small 
compared  to  the  externally  applied  stress  of  120  MPa.  Accordingly,  Og  is  taken 
to  be  zero  in  the  present  models.  We  also  assume  that  the  properties  js  and  xj/ 
are  the  same  at  all  the  grain  boundaries  of  the  single  phase  gamma,  dual 
phase  equiaxed  {y+OL2),  and  fully  lamellar  (Y+a2)  TiAl  models.  Reference 


time,  tE  is  calculated  for  the  fully  lamellar  alloy  and  used  to  normalize  the 
time  for  all  models.  In  the  present  work,  grain  boundary  diffusivity  Db  is 
assumed  to  have  the  same  value  for  the  single  phase  models  and  the  y/y 
grain  boundaries  in  the  dual  phase  equiaxed  models.  Since,  Db  of  the 
gamma  phase  is  higher  than  that  for  the  a2  phase  [133,  134],  we  have  assumed 
that  the  Db  of  the  fully  lamellar  alloy  is  approximately  equal  to  that  for  the 
gamma  phase.  Grain  boundary  properties  of  gamma  TiAl  are  given  in  Table 
6-3.  Although  the  above  assumptions  could  result  in  significant  errors  in 
predicting  absolute  rupture  time,  they  appear  to  be  sufficient  for  an  accurate 
comparison  of  the  susceptibility  of  different  phase  morphologies  to  grain 
boundary  cavitation. 


Table  6-3.  Grain  boundary  properties 


Property 

Value 

Db 

7.5825  X  10-7m2/s  [133] 

Sb 

0.001  |im  [133] 

V 

70°  [131] 

Q. 

4.62x10-30  m3  [133] 

6.3.  Results  and  Discussion 

6.3.1.  Creep  deformation 

Normal  stress  (a^)  contours  for  the  homogeneous  model  are  plotted  in 
Figure  6-4.  It  can  be  seen  that  the  stress  distribution  is  identical  from  grain  to 
grain  in  the  unit  cell.  As  expected,  regions  near  grain  boundary  triple  points 
correspond  to  stress  concentrations.  It  had  been  shown  analytically  that  grain 
boundary  triple  points  have  a  singularity  in  stress  when  grain  boundary 
sliding  is  accommodated  by  power  law  creep  of  the  grain  matrix.  In  the 
present  analysis,  the  maximum  stress  value  at  the  grain  boundary  triple  point 
is  of  the  order  of  two  times  the  externally  applied  stress.  However,  this  error 
will  have  an  insignificant  effect  in  the  results  of  the  present  analysis  since  the 
singularity  is  confined  to  a  very  small  volume. 

Figure  6-5  shows  a  plot  of  normal  stress  contour  in  the  dual  phase  model  in 
absence  of  grain  boundary  sliding  (GBS).  It  is  evident  from  this  figure  that  the 
central  grain,  which  is  of  a2-Ti3Al  phase,  supports  a  low  level  of  stress 
compared  to  that  for  the  adjacent  y  TiAl  grains.  Grain  boundary  triple  points 
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Figure  6-4.  Axial  stress  (Ozz)  distribution  in  single  phase  yTiAl  in  the  presence  of  grain 
boundary  sliding.  The  fringe  levels  defined  on  the  right  side  are  in  units  of  MPa. 
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Figure  6-5.  Axial  stress  (Ozz)  distribution  in  dual  phase  a2+J  TiAl  without  grain 
boundary  sliding.  The  fringe  levels  defined  on  the  right  side  are  in  units  of  MPa. 
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J  and  K  are  found  to  be  in  regions  of  high  stress  concentrations.  However,  the 
four  other  triple  points  (R,  S,  Q  and  P)  in  the  unit  cell  are  relatively  free  from 
stress  concentrations.  Since  the  y  phase  has  a  greater  creep  strength  that  the 
a2  phase  under  present  conditions,  load  is  shed  from  the  grain  to  adjacent 
y  grains  which  results  in  the  low  level  of  at  the  transverse  grain  (a2/  y) 
boundaries.  By  contrast,  transverse  grain  boundaries  between  two  y/y  grains 
(GJ  and  KB)  are  subjected  to  a  high  level  of  02- 

The  features  of  the  Oz  fringe  pattern  for  the  dual  phase  model  in  presence  of 
GBS  are  similar  to  that  of  Figure  6-4,  except  that  the  stress  distribution  inside 
the  central  grain  (a2  phase)  is  different  from  that  in  adjoining  grains  (y 
phase).  Moreover,  grain  boundary  triple  points  J  and  K  have  a  higher  level  of 
stress  concentration  compared  to  other  triple  points  R,  S,  Q  and  P. 

The  normal  stress  fringe  pattern  for  the  FL-1  model  without  GBS  is  illustrated 
in  Figure  6-6.  It  is  interesting  to  note  the  stress  distribution  in  the  a2  laths  in 
this  figure.  Those  a2  lamellae  which  are  inclined  at  an  angle  close  to 
transverse  direction  are  found  to  have  a  Oz  level  almost  identical  to  the 
adjacent  y  laths.  As  the  angle  of  inclination  increases,  stress  in  02  laths  is 
redistributed  and,  as  a  result,  the  adjacent  y  laths  are  under  a  higher  stress 
level.  A  similar  stress  pattern  was  obtained  in  FL-2  model  without  GBS. 

The  O2  distribution  in  the  FL-1  model  with  grain  boundary  sliding  is 
illustrated  in  Figure  6-7.  The  stress  distribution  is  nonuniform  throughout 
the  structure.  It  is  not  evident  from  this  figure  that  the  angular  orientation  of 
the  y  and  a2  phase  boundaries  have  an  effort  on  the  stress  distribution. 
However,  abrupt  changes  in  stress  contours  inside  the  grains  result  from  the 
alternate  change  of  phases.  As  in  the  case  of  the  single  phase  or  dual  phase 
equiaxed  models,  the  transverse  grain  boundaries  are  observed  to  be  the  most 
heavily  stressed  regions.  A  similar  stress  pattern  was  obtained  for  the  FL-2 
model. 

The  normal  stress  distribution  in  FL-1  model  having  y/y  lamellar  boundary 
sliding  is  illustrated  in  Figure  6-8.  Grain  boundary  triple  points  P,  Q,  K  and  B 
are  found  to  be  at  regions  of  stress  concentration.  It  can  also  be  seen  in  Figure 
6-8  that  a  few  regions  along  grain  boundaries  RJ,  JP,  PQ,  QD  and  QK  have 
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Figure  6-6.  Axial  stress  (Gzz)  distribution  in  fully  lamellar  TiAl  (FL-1)  without  grain 
boundary  sliding.  The  fringe  levels  defined  on  the  right  side  of  this  figure  are  in  units  in 
MPa. 
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Figure  6-7.  Axial  stress  (Ozz)  distribution  in  fully  lamellar  TiAl  (FL-1)  in  the  presenee  of 
grain  boundary  sliding.  The  fringe  levels  defined  on  the  right  side  of  this  figure  are  in 
units  in  MPa. 


stress  concentrations  which  resulted  from  the  slide  lines  (y/y  phase 
boundaries)  that  are  constrained  at  these  grain  boundaries.  It  appears  from 
this  figure  that  the  stress  distribution  in  each  grain  is  affected  by  the 
orientation  of  the  laths.  Grains  1  and  6  have  a  relatively  uniform  Oz 
distribution  which  results  from  the  near  transverse  lath  orientation  in  these 
two  grains. 

For  the  normal  stress  distribution  in  FL-2  model  having  y/y  lamellar 
boundary  sliding,  fewer  stress  concentrations  result  for  this  model  compared 
to  that  for  the  FL-1  model.  It  appears  that  the  asymmetric  orientation  of  laths 
allows  for  greater  stress  relaxation  where  the  sliding  phase  boundaries  are 
constrained,  particularly  along  grain  boundary  PQ. 

A  comparison  of  overall  strain  rate  as  a  function  of  strain  for  the  different 
models  is  shown  in  Figure  6-9.  It  is  evident  from  this  figure  that  the  single 
phase  model  without  GBS  exhibits  the  lowest  minimum  strain  rate.  On  the 
other  hand,  FL-1  with  grain  boundary  sliding  exhibits  the  highest  minimum 
creep  rate.  We  note  from  this  figure  that  the  effect  of  grain  boundary  sliding 
is  very  similar  in  nature  from  one  microstructure  to  another.  It  is  also 
evident  from  this  figure  that  the  models  exhibit  a  strain  softening  behavior. 
This  softening  behavior  is  a  result  of  the  8%  increase  in  the  stress  throughout 
the  deformation  process. 

Figure  6-10  shows  the  minimum  strain  rates  of  FL-1  and  FL-2  models  with 
y/y  lath  boundary  sliding.  It  can  be  seen  in  this  figure  that  the  minimum 
strain  rate  is  higher  for  the  FL-1  model  by  a  factor  of  4  as  compared  to  that  in 
FL-2  model.  A  separate  calculation  based  on  a  shear  crack  model  and 
subsequently  modified  for  a  random  microstructure  was  performed  as  an 
alternate  method  to  estimate  the  effect  of  y/y  phase  boundary  sliding  in  these 
two  models.  In  spite  of  some  sources  of  error  in  applying  the  shear  crack 
model  to  phase  boundary  sliding,  we  find  that  the  trend  in  the  results 
obtained  using  a  shear  crack  model  is  similar  to  that  obtained  in  present 
modeling  using  FEM.  We  note  that  this  shear  crack  approximation  is  good 
for  a  rough  estimation  of  phase  boundary  sliding  but  a  detailed  finite  element 
analysis  must  be  carried  out  in  order  to  gain  more  insight  into  the  problem. 
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Figure  6-8.  Axial  stress  (Gzz)  distribution  in  fully  lamellar  TiAl  (FL-1)  in  the  presence  of 
y/y  phase  boundary  sliding.  The  fringe  levels  defined  on  the  right  side  of  this  figure  are  in 
units  in  MPa. 
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Figure  6-9.  Comparison  of  overall  strain  rate  vs.  strain  for  different  TiAl  microstructures. 
Models  which  include  grain  boundary  sliding  are  indicated  by  w/GBS. 
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Figure  6-10.  Comparison  of  overall  strain  rate  vs.  strain  in 
FL  TiAl  for  lath  boundary  sliding  in  the  absence  of  grain 
boundary  sliding. 
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Facet  stress  enhancement  factors  were  calculated  for  all  the  models  and  are 
compared  in  Figure  6-11.  Anderson  and  Rice  used  a  self  consistent  approach 
to  study  the  effect  of  constraints  on  diffusive  creep  cavity  growth  for  f.c.c. 
metals.  They  found  a  value  of  facet  stress  enhancement  factor,  f^  =  1.67  for  n 
=  5  for  3D  Wagner-Seitz  cells  when  hexagonal  facets  are  normal  to  externally 
applied  stress.  Rice  has  shown  from  2D  analysis  that  fa  for  a  homogeneous 
single  material  with  freely  sliding  grain  boundaries  is  approximately  1.5.  We 
obtained  f^  =  1.49  for  the  single  phase  model,  as  shown  in  Figure  6-11. 

The  single  phase  model  exhibits  a  greater  value  for  f^y,  than  that  obtained  at 
the  y/y  grain  boundaries  of  dual  phase  equiaxed  model.  However,  f^  at  the 
a2/y  grain  boundaries  of  the  dual  phase  model  is  greater  than  the  facet  stress 
enhancement  factor  of  single  phase  model.  Grain  boundary  sliding  in  the  FL- 
1  and  FL-2  models  yields  a  lower  value  of  f^y  when  compared  to  single  phase 
and  dual  phase  equiaxed  models.  We  observe  that  the  minimum  strain  rate 
is  enhanced  due  to  grain  boundary  sliding  in  the  single  phase  and  dual  phase 
equiaxed  structures  by  factors  of  2.8  and  3.3  respectively.  On  average,  the 
stress  in  the  y  grains  is  higher  than  that  supported  by  the  a2  grains.  In  the 
single  phase  structure,  the  mean  stress  in  the  y  grains  is  lower  resulting  in  a 
minimum  strain  rate  as  compared  to  the  dual  equiaxed  model  by  a  factor  of 
1.6.  This  observation  is  consistent  with  the  experimental  findings  that  under 
similar  creep  conditions,  a  duplex  TiAl  alloy  indeed  exhibits  a  higher  creep 
strain  rate  than  that  for  the  corresponding  single  phase  y  microstructure  by  a 
factor  of  1.83. 

Of  all  of  the  present  forms  of  TiAl  considered  in  a  present  work,  the  FL 
microstructure  was  found  experimentally  to  possess  the  highest  creep 
strength  and  fracture  toughness.  This  finding  can  only  be  consistent  with  the 
present  model  results  if  grain  boundary  sliding  is  inhibited  in  the  fully 
lamellar  microstructure  while  relatively  easy  in  the  equiaxed  microstructures. 
Two  different  types  of  lamellar  grain  boundaries  had  been  found  in  FL  TiAl 
alloy.  One  with  a  composition  of  Ti-  44%,  is  found  to  have  serrated  grain 
boundaries.  It  was  observed  that  the  alloy  having  serrated  grain  boundaries 
exhibits  a  lower  creep  strain  rate  compared  to  the  alloy  with  straight  grain 
boundaries  by  a  factor  of  8  and  concluded  that  the  serrated  nature  of  the  grain 
boundary  is  the  primary  reason  for  higher  creep  strength.  This  is  in  accord 


with  the  present  FE  results  given  that  the  straight  grain  boundaries  undergo 
sliding  during  creep  but  the  serrated  grain  boundaries  are  generally  prevented 
from  sliding.  A  comparison  of  Figures  6-9  and  6-10  indicates  that  the  FL 
alloys  would  exhibit  an  overall  strain  rate  in  presence  of  y/ y  lath  boundary 
sliding  that  is  about  two  orders  of  magnitude  greater  than  that  for  the 
equiaxed  microstructures  in  the  presence  of  grain  boundary  sliding.  This 
result  appears  unrealistic  as  the  overall  strain  rates  obtained  from  the  lath 
boundary  sliding  models  are  more  than  hundred  times  greater  than  the 
typical  minimum  creep  strain  rates  observed  experimentally  for  LF  TiAl 
alloys  under  similar  temperature  and  stress  conditions.  Moreover,  it  was 
observed  experimentally  that  the  density  of  superdislocations  <101]  in  y 
lamella  do  not  change  during  the  course  of  creep  deformation.  Therefore,  it 
was  concluded  that  the  superdislocations  do  not  favor  any  extra  deformation 
mode  along  the  y/y  laths.  Thus  it  is  doubtful  that  many  superdislocations  are 
able  to  reach  the  lath  boundaries  to  dissociate  under  the  creep  conditions 
imposed.  Based  on  this  data  and  our  model  results,  it  is  unlikely  that  y/ y 
phase  boundary  sliding  occurs  to  any  great  degree  during  the  creep 
deformation  of  lamellar  TiAl  alloys. 

6.3.2.  Creep  constrained  cavitation 

6.3.2.I.  Cavitation  in  fully  lamellar  and  dual  phase  equiaxed  models 
Figure  6-12  shows  the  normalized  axial  stress  distribution  in  the  fully 
lamellar  model  (FL-1)  in  the  absence  of  GBS.  It  is  observed  that  az  laths  that 
are  more  closely  aligned  with  the  applied  stress  direction  redistribute  more 
stress  to  the  adjacent  gamma  laths.  This  results  in  inhomogeneous  stress 
distributions  along  grain  boundaries  RS  and  PQ.  However,  the  stress  level 
discontinuities  are  more  extreme  at  boundary  RS  than  those  at  PQ.  This  is 
apparently  due  to  the  more  isostress  orientation  of  the  lamellae  in  grain  5  as 
well  as  the  proximity  of  RS  to  the  cavitating  facet. 

The  distribution  of  normalized  axial  stress  in  FL-1  model  in  the  presence  of 
GBS  is  illustrated  in  Figure  6-13.  It  is  evident  from  this  figure  that  the 
transverse  grain  boundaries  RS  and  KB  are  heavily  stressed  regions  in  the 
microstructure.  Cavitation  at  grain  boundary  OH  again  redistributes  stress 
from  grain  1  to  adjacent  grains.  However,  a  high  level  of  stress  at  grain 
boundary  RS  ensues  as  a  result  of  the  load  shedding  of  the  az  laths  at  this 
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Figure  6-12.  Normalized  axial  stress  (Gz/^)  distribution  in  the  fully  lamellar  TiAl  (FL-1) 
during  creep  cavitation  at  central  grain  boundary  OH  without  GBS.  It  should  be  noted  that 
the  fringe  levels  in  this  figure  are  different  from  those  in  previous  figures. 
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Figure  6-13.  Normalized  axial  stress  {<5jL)  distribution  in  the  fully  lamellar  TiAl  (FL-1) 
during  creep  cavitation  at  central  grain  boundary  OH  in  the  presence  of  GBS. 


grain  boundary  and  due  to  the  sliding  of  the  inclined  grain  boundaries.  It  is 
interesting  to  note  that  the  transverse  grain  boundary  KB  is  stressed  to  a 
higher  level  than  either  boundary  GJ  or  PQ  although  these  boundaries  have 
lath  orientations  that  are  less  characteristic  of  an  isostress  arrangement.  It 
appears  that  the  higher  stress  at  boundary  KB  is  due  to  the  fact  that  it  lies  in 
the  path  of  the  highest  stresses  which  arise  from  the  growth  of  the  cavities  at 
grain  boundary  OH. 

The  development  of  creep 
cavitation  damage  co=a/b  as  a 
function  of  normalized  time 
t/tg  is  illustrated  in  Figure  6- 
14,  for  the  present  models.  It 
is  evident  from  this  figure 
that  grain  boundary  sliding 
substantially  reduces  the 
cavity  coalescence  time  for 
all  the  microstructures 
considered. 

The  two  fully  lamellar 
models,  FL-1  and  FL-2 

yielded  identical  results  for  cavitation  in  the  presence  of  GBS.  This  similarity 
indicates  that  the  effect  of  lamellae  orientation  on  creep  cavitation  is  not  as 
pronounced  in  the  presence  of  grain  boundary  sliding.  In  the  absence  of  GBS, 
the  FL-2  model  is  found  to  be  stronger  than  the  FL-1  model.  The  important 
point  here  is  that  a  difference  in  cavity  coalescence  time  for  the  two  FL 
models  is  only  evident  when  GBS  is  suppressed.  This  is  likely  the  case  for  the 
FL  microstructure  which  normally  exhibits  serrated  grain  boundaries.  The 
present  results  indicate  that  cavities  at  a  grain  boundary  with  intersecting  lath 
orientations  which  are  isostress  in  nature  will  tend  to  grow  slower  compared 
to  the  case  where  the  intersecting  lath  orientations  are  closer  to  an  isostrain 
orientation.  Lath  orientations  that  are  nearly  parallel  to  the  applied  stress 
direction  in  the  zone  adjacent  to  the  cavitating  grain  boundary  result  in  load 


t/t^ 


Figure  6-14.  Plot  of  creep  damage  at  the  facet  #1  in 
FL  and  DX  models  versus  normalized  time. 


shedding  from  OLi  laths  to  y  laths  which  then  produces  less  constraint  on  the 
cavitation  process. 

6. 3. 2. 2.  'Effect  of  lamella  thickness 

Effect  of  lath  width  on  creep  constrained  grain  boundary  cavitation  in  the 
fully  lamellar  TiAl  intermetallics  was  also  investigated.  Two  models  FL-lb 
and  FL-2b  having  l/d=0.058  are  analyzed  here.  1  is  the  lath  thickness  and  d  is 
the  grain  size  with  d=4Ro.  The  FL-1  and  FL-2  models  have  1/ d=0.029.  The 
meshes  used  for  FL-lb  and  FL-2b  models  are  identical  with  those  used  for  FL-1 
and  FL-2  models  except  that  the  laths  have  twice  the  thickness. 

Figure  6-15  shows  the 
comparative  growth  of  creep 
damage  in  the  fully  lamellar 
models  FL-1,  FL-2,  FL-lb  and 
FL-2b.  Cavity  growth  as  a 
function  of  normalized  time 
t/tE  is  calculated  at  facet  OH 
in  models  FL-1  and  FL-lb 
and  at  facet  EF  for  FL-2  and 
FL-2b  models.  It  is  evident 
from  this  figure  that  the 
development  of  cavitation 
damage  for  all  the  models  is 
practically  identical  in  the 

presence  of  the  GBS.  This  is  due  to  the  fact  that  when  grain  boundary  sliding 
takes  place,  the  regions  around  the  grain  boundaries  carry  the  bulk  of  the  load 
and  deform  heavily  which  overrides  the  effect  the  laths  have  on  stress 
redistribution  in  the  adjoining  grains.  In  the  absence  of  grain  boundary 
sliding,  the  thicker  lath  width  resulted  in  faster  cavity  growth  rates  for  both 
the  FL  models.  However,  these  differences  in  cavity  growth  rates  are  small 
and  for  all  practical  purposes  may  be  deemed  negligible.  It  is  also  evident 
from  this  figure  that  the  cavity  growth  rate  is  consistently  dependent  on  the 
lamella  orientation  in  the  absence  of  grain  boundary  sliding.  Laths  in  the 
regions  immediately  adjacent  to  the  cavitating  facets  in  the  FL-2  and  FL-2b 


t/te 

Figure  6-15.  Plot  of  creep  damage  at  the  facet  #1  in 
FL-1  and  FL-2  models  versus  normalized  time. 
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models  are  of  isostress  orientation  providing  more  constraint  on  cavitation  as 
compared  to  that  for  the  FL-1  and  FL-lb  models. 

It  appears  from  experimental  observations  that  a  thicker  lath  width  can  result 
in  a  lower  degree  of  lamella  interlocking  at  the  grain  boundaries  [5]. 
Therefore,  grain  boundary  sliding  mobility  could  also  be  higher  in  a 
microstructure  having  thicker  lath  width.  It  follows  that,  the  effect  of  thicker 
lath  width  on  grain  boundary  cavitation  could  be  detrimental  due  to  more 
stress  redistribution  from  sliding  grain  boundaries.  A  good  understanding  of 
the  relationship  between  lath  thickness  and  grain  boundary  sliding  mobility 
in  TiAl  requires  further  study. 


w/GBS 

Facet  #1  (FL-iy 
Facet  #2  (FL-1) 
Facet  #1  (FL-2) 
Facet  #2  (FL-2) 

Facet  #1  (FL-1) 
Facet  #2  (FL-1) 
Facet  #1  (FL-2) 
Facet  #2  (FL-2) 


6. 3.2. 3.  Effect  of  interaction  between  cavitatin^  facets 

Figure  6-16  shows  the 

development  of  creep 

damage,  omega  due  to  the 

interaction  between  the 

cavitating  facets  1  and  2  as  a 

function  of  normalized  time, 

t/tE  for  the  fully  lamellar 

models,  FL-1  and  FL-2.  As  in 

the  previous  figure,  GBS 

results  in  a  significant 

reduction  in  rupture  time  for 

the  present  FL  models.  It  is  ^ 

Figure  6-16.  Plot  of  creep  damage  at  the  facets  #1 
evident  here  that  the  growth  #2  in  FL-1  and  FL-2  models  versus  normalized 

of  cavities  at  the  two  facets  are  time, 
nearly  identical  in  all  the 

models.  Cavitation  damage  in  one  facet  redistributes  stresses  to  its  adjacent 
transverse  grain  boundaries.  Thus,  a  positive  interaction  between  the 
cavitating  facets  1  and  2  takes  place  which  results  in  nearly  identical  cavity 
growth  in  both  the  facets.  As  in  the  previous  section,  the  FL-2  model  is  found 
to  be  more  resistant  to  cavity  growth  than  the  FL-1  model  in  the  absence  of 
GBS.  The  difference  in  rupture  time  for  all  the  models  with  GBS  is  not  very 
large  and  hence  can  be  considered  as  identical  for  all  practical  purposes. 
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The  creep  damage  due  to  the 
interaction  between  the 
cavitating  facets  1  and  3  in 
the  FL  models  are  plotted  as 
a  function  of  normalized 
time,  t/tE,  in  Figure  6-17. 

6. 3.2. 4.  Effect  of  cavitating 
facet  density 

Mesh  used  for  studying  the 
effect  of  lower  cavitating 
facet  density  in  FL  TiAl  is 
shown  in  Figure  6-18.  The 
boundary  conditions  used  for  time, 
this  bigger  unit  cell  are 

similar  to  those  applied  to  the  smaller  unit  cell  OACF. 


w/GBS 
Facet  #1  (FL-1) 


- Facet  #3  (FL-1) 

- Facet#!  (FL-2) 

- Facet  #3  (FL-2) 


- Facet  #1  (FL--1) 

- Facet  #3  (FL-1) 

- Facet  #1  (FL-2) 

- Facet  #3  (FL-2) 


% 

Figure  6-17.  Plot  of  creep  damage  at  the  facets  #1 
and  #3  in  FL-1  and  FL-2  models  versus  normalized 


The  effect  of  cavitating  facet  density 
on  creep  cavitation  damage  is 
shown  in  Figure  6-19.  The  FL-1 
model  with  higher  cavitating  facet 
density  yielded  faster  cavity  growth 
for  both  with  and  without  the  GBS. 
The  FL-1  model  with  p  =  0. 03186 
exhibited  rapid  cavitation  damage 
growth  in  the  absence  of  GBS  from 
the  very  onset  of  cavitation. 
However,  the  two  FL-1  models 
with  p=0. 01416  and  p  =  0. 03186 
yielded  identical  damage  growth 
until  the  cavities  reach  a  critical 
stage  at  (o=0.42  in  the  presence  of 
GBS.  From  this  time  onward,  the 
FL-1  model  with  p=0.03186  yielded 


Figure  6-18.  FE  mesh  for  the  FL-1  and  FL- 
2  models  having  r=0.01416.  The  a2  laths  are 
indicated  by  darker  shading. 


much  faster  damage  growth  till  the  cavities  coalesce.  It  should  be  noted  here 
that  an  increase  of  p  by  a  factor  of  2.25  only  reduces  the  rupture  time  by  a 


factor  of  1.28  and  1.23  for  the  FL-1  model  with  and  without  the  GBS 
respectively. 


6.4.  Conclusion 

The  following  conclusions  are  drawn  from  the  present  study  of  creep 
deformation  and  cavitation  in  fully  lamellar  and  dual  phase  equiaxed  (a2+y) 
TiAl  microstructures. 


It  is  found  that  inherent 
creep  strength  of  the  fully 
lamellar  TiAl  is  due  to  the 
nature  of  its  serrated  grain 
boundaries  which  inhibits 
boundary  sliding  during 
creep  deformation.  It  is  also 
found  that  the  y/y  boundary 
sliding  does  not  take  place  to 
a  great  extent  during  the 
creep  deformation.  These 
findings  are  consistent  with 
the  experimental 

observations  published  in 
the  literature. 


Figure  6-19.  Effect  of  cavitating  facet  density  on 
creep  damage  of  the  fully  lamellar  models  versus 
normalized  time. 


The  present  model  results  indicate  that  the  dual  phase  equiaxed 
microstructure  exhibits  the  greatest  resistance  against  cavitation  when  cavity 
growth  occurs  at  the  y/y  grain  boundary  adjacent  to  the  02  grain.  This  occurs 
in  the  presence  of  grain  boundary  sliding.  Conversely,  this  microstructure  is 
the  most  susceptible  to  cavitation  damage  when  cavitation  occurs  at  an 
inclined  angle  to  a  neighboring  a2  grain.  The  effect  of  lamellae  orientation  on 
creep  cavitation  damage  in  the  fully  lamellar  microstructures  is  found  to  be 
insignificant  in  the  presence  of  grain  boundary  sliding.  However,  isostress 
orientation  of  laths  in  the  grains  adjacent  to  the  cavitating  grain  boundary  is 
found  to  delay  the  cavity  growth  and  thereby  increases  the  rupture  life. 
However,  in  real  microstructures  with  random  lamellae  orientation,  this 


beneficial  effect  may  not  be  seen.  It  has  been  observed  in  the  present  study 
that  the  lamellae  width  does  not  affect  the  cavity  growth  process  as  long  as  the 
grain  boundary  mobility  is  the  same  in  all  the  different  fully  lamellar 
microstructures.  The  effect  of  interaction  between  the  cavitating  facets  in  the 
fully  lamellar  models  is  generally  positive  when  cavitating  boundaries  are  on 
neighboring  grains.  However,  a  negative  interaction  effect  can  take  place  for 
certain  cavitating  facet  configuration. 

Based  on  the  present  results,  the  relatively  long  rupture  times  observed 
experimentally  for  fully  lamellar  TiAl  may  be  primarily  attributed  to  the 
suppression  of  grain  boundary  sliding  in  this  microstructure. 
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3.  "Investigation  on  the  Creep  Behavior  of  Spray  Formed  y-TiAl  Alloys"  B. 
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Composite  using  Finite  Element  Method  and  X-Ray  Diffraction"  S.  Ho 
and  E.J.  Lavernia,  presented  at  the  TMS/ASM  International  Symposium 
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Lavernia,  presented  at  the  4th  International  Conference  on  Powder 
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"Characterization  of  Low-Cycle  Fatigue  Damage  on  Inconel  718  by  Laser 
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"Microstructure  and  Fracture  Behavior  of  a  Spray  Atomized  and 
Deposited  Nickel  Aluminide  Intermetallic",  D.  Lawrynowicz,  X.  Liang, 
T.S.  Srivatsan  and  E.J.  Lavernia,  presented  at  the  Second  International 


Symposium  on  Fatigue  and  Fracture  of  Ordered  Intermetallic  Materials. 

October  2-6,  1994,  Rosemont,  Illinois. 

14.  "Analysis  of  a  Dendritic  NisAl  Alloy  by  X-Ray  Diffraction",  K.  J.  C.  Chou 
and  J.  C.  Earthman,  presented  at  the  International  Conference  on  High- 
Temperature  Intermetallics,  May  16-19,  1994,  San  Diego,  California. 


Al.1.2.  Posters 

1.  "Finite  Element  Analysis  of  Primary  Creep  in  Fully  Lamellar  yTiAl 
Intermetallics",  A.  Chakraborty  and  J.C.  Earthman,  postered  at  the  MRS 
Spring  Meeting.  April  8-12,  1996,  San  Francisco,  California. 
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1.  Interactions  with  Dr.  Lee  Semiatin  (Wright  Laboratories,  WL/MLLN 
Wright  Patterson  AFB,  OHIO  45433).  During  the  AFOSR  funded  period,  we 
actively  exchanged  results,  including  both  processing  details  as  well  as 
compressive  creep  data  obtained  in  our  AFOSR  funded  research  program 
with  Dr.  Semiatin.  Dr.  Semiatin's  comments  on  our  data  were  very  useful  in 
the  interpretation  of  the  results  obtained  in  the  present  AFOSR  program. 

2.  Interactions  with  Dr.  Don  Larsen  (Howmet  Corporation,  Operhall 
Research  Center,  1500  South  Warner  Street,  Whitehall,  MI  49461-1895).  Dr. 
Larsen  has  been  actively  working  with  cast  y-TiAl  for  a  number  of  years. 
Presently,  the  Howmet  Corporation  is  actively  working  with  a  variety  of 
compositions  for  elevated  temperature  aerospace  applications.  Moreover,  the 
Howmet  Corporation  is  actively  developing  spray  forming  technology  as  a 
manufacturing  approach  for  the  spraying  of  superalloy  rings.  Since  the  pilot 
plant  spray  facility  that  is  available  at  Howmet  is  not  designed  for  the  spraying 
of  Ti,  Dr.  Larsen  expressed  a  strong  interest  in  our  AFOSR  funded  program. 
During  the  AFOSR  funded  period,  we  successfully  spray  formed  one  of  the 
titanium  aluminide  compositions  provided  by  Howmet  Corporation.  The 
spray  formed  titanium  aluminide  was  then  sent  to  Howmet  for 
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3.  Interactions  with  Dr.  Shyh-Chin  Huang  (General  Electric  Research  and 
Development  Center,  Building  K-1,  Room  MB  275,  P.O.  Box  8,  Schenectady, 
NY  12301).  Dr.  Huang  has  extensive  experiences  in  the  study  of  y-TiAl  alloys. 
He  has  been  well  known  for  the  development  of  Ti-48Al-2Nb-2Cr  alloys.  Our 
interactions  with  Dr.  Huang  included  successful  spray  forming  of  Ti-42A1- 
2Nb-2Cr  alloys  supplied  by  General  Electric  Research  and  Development 
Center,  characterization  of  spray  formed  Ti-42Al-2Nb-2Cr  in  General  Electric, 
and  exchange  of  the  preliminary  results. 
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